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INTRODUCTION 


Whereas  transcriptional  gene  re  gulation  is  well  studied,  posttranscr  iptional  gene  control  is  poorly 
understood.  Yet,  em  erging  evidence  indicates  that  posttranscriptiona  1  control  by  RNA  binding  proteins 
(RBPs)  and  m  icroRNAs  (m  iRNAs)  are  im  portant  key  re  gulators  of  gene  expression  of  m  any  cancer 
related  genes.  The  RBP,  HuR,  is  a  m  aster  regulator  of  many  early  response  genes  in  cancer.  Of  the  six 
classical  acquired  traits  first  pr  oposed  by  Hanahan  and  Weinberg  that  m  alignantly  tran  sformed  cells 
develop,  HuR  has  been  demonstrated  to  control  genes  in  multiple  areas. 

The  purpose  of  this  research  is  to  define  the  in  vivo  mRNA  targets  of  the  RNA  binding  protein, 
HuR  using  RNA  im munoprec i p i tat i o n s  applied  to  m  icroarray  chip  s  (  RIP-Chip)  in  estrogen  positiv  e 
(ER+)  and  estrogen  negative  (ER-)  breast  ca  ncer.  Additionally,  we  will  test  the  hypothesis  that  HuR 
increases  cell  proliferation  by  regulatin  g  mRNAs  involved  in  breast  oncogenesis  at  the 
posttranscriptional  level.  We  will  do  this  by  alter  ing  HuR  levels  in  breast  ca  ncer  cell  lines  using 
lentiviral  short  hairpin  RNAi  to  knock  down  and  also  over-express  HuR  and  assay  turn  or  formation  by 
injection  into  athymic  nude  mice  (orthotopic  model). 
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BODY  OF  WORK 


The  funding  from  the  DOD  resulted  in  two  public  ations,  w  hich  are  listed  in  the  Appendix  and 
References:  Calaluce  et  al  and  Gubin  et  al.  In  the  first  publication,  Calaluce  et  al,  we  described  how  we 
developed  methods  of  RIP-Chip  and  applied  them  to  both  ER-  and  ER+  breast  cancers,  MD  A-MB-23 1 
and  MCF-7,  respectively.  In  Figure  1 ,  (Calaluce  et  al)  we  see  that  Hu  R  RIP  is  able  to  identify  a  known 
HuR  target,  [1-actin  mRNA  robustly  in  both  cell  types.  W  e  then  extended  these  RIP  -Chip  techniques  to 
genome  wide  arrays.  In  Figure  2  (Calaluce  et  al),  re -produced  below,  we  see  a  statistical  com  pilation  of 
three  independent  immunoprecipitations  with  matching  controls  in  both  cell  lines.  There  are  three  classes 
of  HuR  ass  ociated  mRNAs:  those  associated  with  HuR  only  in  MD  A-MB-23 1  cells  (upper  left),  those 
only  with  MCF-7  (lower  right)  and  mRNAs  that  are  HuR  associated  with  both  cells  types. 


MDA  HuR  T argets  &  MCF  HuR  Targets  (Log2  scale) 


We  perfonned  Gene  Ontology  (GO)  classification  of  these  genes 
are  suggested  putative  novel  m  embers  to  the  acquired 
capabilities  m  odel  of  cancer  as  proposed  by  Hanahan  and 
Weinberg  (  Figure  3  ,  Calaluce  e  t  al.).  The  o  riginal  RIP-  Chip 
results  were  confir  med  using  RT-  PCR  for  t  wo  target  genes, 
calmodulin  2  (  CALM2)  and  CD9  (Figure  4  ,  Calaluce  et  al). 
These  genes  were  chosen  since  they  have  known  cancer 
associations  and  were  highly  re  presented  as  HuR  t  arget  mRNA 
in  both  cell  types.  As  se  en  in  Figure  4,  HuR 

immunoprecipitation  efficiently  reco  vers  an  enrichm  ent  of  both 
CALM2  and  CD9  as  compared  to  isotype  controls. 

Wever  ifiedt  heRI  P-Chip  result  susi  ng  a  second 
independent  method,  biotin  pull  downs.  In  biotin  pull  downs,  the 
relevant  portion  of  the  target  m  RNA  is  in  vitro  transcribed  with 
biotin-UTP  and  this  extracted  mRNA  is  then  mixed  with  cellular 
lysate  which  contains  HuR  prot  ein.  The  complex  is  recovered 
and  then  analyzed  by  Western.  As  seen  in  the  results  in  Figure  5 
(Calaluce  et  al),  both  CALM2  and  CD9  have  predicted  HuR 
binding  sites  in  the  3  ’  UTR  and  bi  otin  pull  downs  indicate  that 
HuR  interacts  with  only  the  3’  UTR  but  not  any  ot  her  region  of 
the  mRNAs.  W e  conclude  that  by  two  independen  t  methods,  RIP-Chip  and  biotin  pull  downs,  HuR  is 
indeed  physically  associating  with  CALM2  and  CD9  mRNAs. 

Physical  association  between  HuR  a  nd  mRNA  targets  does  not,  however,  prove  functionality.  In 
order  to  test  whether  H  uR  is  contr  oiling  expre  ssion  of  the  se  two  targe  t  mRNAs,  we  m  ade  mutant  cel  1 
lines,  MD  A-MB-23 1  and  MCF-7  which  over-expre  ssed  or  under-expressed  HuR.  As  seen  in  Figures  6 
and  7  (Calaluce  et  al),  using  lentiviral  vectors  whic  h  express  shRNA  HuR,  we  were  able  to  knock-down 
HuR  by  ove  r  90%  in  both  cell  lines.  W  e  also  devel  oped  stably  transfected  clones  which  over-expressed 
epitope  (H  A)  tagged  HuR  by  different  am  ounts  in  both  MD  A-MB-23 1  and  MCF-7  cells.  However, 
unexpectedly,  over-expression  of  Hu  R  in  MB-231  cells  resulted  in  decreases  in  both  CD9  mRNA  and 
protein  ( Figure  6  ,  Calaluce  et  al)  but  the  leve  Is  of  CALM2  were  not  decreased  but  rather  increased. 
However,  the  results  in  MCF-7  cells  were  as  predicted:  HuR  over-expression  resulted  in  increased  in  both 
CD9  and  CALM2  m  RNA  and  protein  (  Figure  7,  Calaluce  et  al).  Therefore,  we  conclude  that  altering 
HuR  levels  affect  the  same  genes  differently,  depending  upon  cellular  milieu. 

We  believe  that  th  e  importance  of  this  work  lie  s  in  two  are  as:  firstly,  we  were  able  to  identif  y 
cancer-associated  genes  as  novel  HuR  targets,  which  was  previously  unknown.  Secondly,  we  believe  that 
the  RIP-Chip  m  ethods  m  ay  have  uncovered  novel  can  cer  genes.  There  were  m  any  genes  which  were 
strongly  recruited  to  the  HuR  pelle  t  in  both  cell  types,  which  do  not  have  any  known  cancer  functions. 
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Figure  2:  HuR  puli  downs  reveal  unique 
sets  of  mRNA  targets  in  breast  cancer  cell 
lines.  HuR  pull  downs  were  performed  from 
MCF-7  or  MDA  MD  231  lysates  using  HuR 
antibody  (3A2)  and  IgGI  isotype  control 
and  hybridized  to  lllumina  Sentrix  arrays 
(47,000  genes).  Control  signals  are 
subtracted.  Above  results  represent 
cumulative  data  from  12  different  arrays. 
Experiments  were  done  n=3.  Scales  are 
log2.(  From  Calaluce  et  al.) 


6 


Since  HuR  generally  af  fects  genes  by  increasing  th  eir  expression,  these  m  ay  well  represent  genes  which 
are  recruited  to  the  polysomes  to  make  more  protein. 

To  more  fully  investigate  the  functional  conseq  uences  of  altering  HuR  expression  in  ER-  breast 
cancer,  we  injected  MD  A-MB-231  cells  with  increas  ed  levels  of  HA-HuR  expression  in  athym  ic  nude 
mice.  The  results  are  described  in  th  e  second  publication,  G  ubin  et  al.  In  Figure  1  (Gubin  et  al),  we  see 
that  MB-231  cells  which  over-expre  ss  HA-HuR  grow  significantly  faster  in  vitro  with  alterations  in  cell 
cycle  kinetics.  Unexpectedly,  when  they  were  introduced  orthotopically  into  mammary  fat  pads  of  female 
athymic  nude  m  ice,  they  produced  m  uch  s  mailer  turn  ors  as  com  pared  to  em  pty  vector  (EV)  controls 
(Figure  ID  ,  Gubin  et  al).  This  surprisingl  y  r  esult  was  r  eplicated  with  a  second  experim  ent,  which 
included  parental  MDA-MB-23 1  cells  (  Figure  2  A,  Gubin  et  al).  There  is  a  90%  reduction  both  in  turn  or 
size  and  mass  of  HA-  HuR  turn  ors.  Turn  ors  dissected  from  ani  mals  contained  live  cells  w  hich  are 
undifferentiated  adenocarcinoma  (Figure  2D).  We  performed  gene  profiling  and  identified  genes  which 
were  up-regulated  in  HA-HuR  tumors  (Table  1,  Gubin  et  al.).  There  were  m  any  genes  of  interest  and  we 
focused  on  known  HuR  target  genes,  such  as  thro  mbospondinl  (TSP1),  and  VEGF.  Additionally  ,  TSP2 
was  also  identified,  though  this  has  not  bee  appreciated  to  be  a  HuR  target.  W  e  asked  whether  there  are 
alterations  of  TSP1  and  VEGF  expr  ession  in  HA-HuR  tumors.  There  are  increases  in  TSP1  expression  at 
both  the  mRNA  and  protein  levels,  but  decreases  in  VEGF  ( Figure  3,  Gubin  et  al.).  We  verified  that  both 
TSP1  and  VEGF  inR  NA  interact  with  HuR  in  HA-  HuR  turn  ors  and  dem  onstrated  that  HuR  over¬ 
expression  results  in  in  creased  stabilization  of  TSP1  but  not  VEGF  mRNA  half-life  (  Figure  4,  Gubin  et 
al.).  Th  ese  results  are  significan  t  because  TSP1  is  know  n  to  be  a  potent  inhib  itor  of  angiogenesis, 
whereas,  VEGF,  is  a  potent  pro-angiogenic  gene.  W  e  ruled  out  putative  explanati  ons,  such  as  alterations 
in  senescence,  proliferation  or  a  poptosis.  The  putative  m  echanism  seems  to  be  a  drastic  reduction  (75%) 
in  blood  vessels  in  HA-HuR  turn  ors.  We  conclude  that  HuR  over-expression  does  indeed  result  in  m  uch 
smaller  tumors  which  has  significant  reductions  in  blood  vessels. 

Finally,  we  have  used  lentiviral  vectors  to  decrease  HuR  levels  in  MDA-MB-23 1  cells  by  making 
stably  infected  clones  (s  ee  Preliminary  Data  Figure  1  ).  When  these  c  ells  were  injected  into  athym  ic 
nude  m  ice,  there  were  growth  reductions  of  95%  in  turn  ors  with  lower  HuR  levels,  as  compared  to 
lentiviral  vector  con  trols.  These  surprising  results  are  c  onfusing  and  seem  ingly  contradict  our  published 
results  in  Gubin  et  al.  However,  there  is  confusion  in  the  literature  whether  HuR  is  a  tumor  suppressor  or 
oncogene  (see  below). 


Preliminary  Data  Figure  1:  MDA-MB-231  cells  infected 
with  a  lentivirus  expressing  a  shRNA  knocking  down  HuR 
(LL  HuR  shRNA)  showed  significantly  reduced  tumor 
volume  (mm3)  and  growth  starting  at  seven  weeks  post¬ 
inoculation  and  continuing  for  fourteen  weeks  when 
compared  to  MDA-MB-231  infected  with  control  lentivirus 
(LL  control).  Five  animals  per  group  were  used;  p<0.05. 


7 


RECOMMENED  CHANGES  FOR  FUTURE  WORK 


Our  findings  for  this  application  m  ust  be  tem  pered  by  two  salient  facts:  (1)  athymic  nude  m  ice 
lack  a  normal  imm  une  system,  which  human  patients  po  ssess;  (2)  we  worked  with  only  one  type  of  ER- 
breast  cancer,  MDA-MB-23 1 .  Therefore,  presently,  it  is  not  clear  whethe  r  our  findings  can  be 
extrapolated  to  other  breast  cancer  types.  The  initial  confusion  in  resu  Its,  suggesting  that  over-expression 
and  under-expression  of  HuR  can  both  reduce  tumor  growth  in  animal  models,  is  reflected  in  the  literature 
over  whether  HuR  is  a  turn  or  suppressor  or  oncogene  s.  Heinonen  and  colleagues  demonstrated  a  direct 
correlation  between  increased  cytoplasm  ic  HuR  and  worse  prognosis  in  breast  cancer  (Heinonen  et  ah). 
Ortega  et  al.  arrived  at  the  opposite  conclusion  (Ortega  et  al).  Therefore,  we  conclude  that  further  study  is 
warranted  to  better  understand  the  role  of  HuR  in  breast  cancer. 

For  future  work,  we  would  s  uggest  two  broad  appro  aches.  First,  we  would  suggest  using  a 
spontaneously  occurring  m  ouse  breast  cancer  model,  su  ch  as  4T1  Balb/c  turn  ors  in  syngeneic  anim  als. 
Syngeneic  anim  als  have  norm  al  immune  system  s.  One  could  m  ake  4T 1  clone  s  which  over-express  or 
under-express  HuR  and  inject  orthotopi  cally  into  Balb/c  mice  and  measure  tumor  growth  and  metastasis. 
Secondly,  we  have  made  a  condi  tional  HuR  knock-out  mouse.  This  mouse  could  be  crossed  with  a 
genetically  engineered  mouse  m  odel  of  breast  cancer,  for  exam  pie,  MMTV-PyMT.  Such  a  m  odel  would 
enable  us  to  genetically  ablate  HuR  and  ask  the  quest  ions  whether  this  in  terferes  with  bre  ast  canc  er 
growth  and/or  metastasis. 

PERSONAL  WHO  RECEIVED  SUPPORT 

Dr.  Ulus  Atasoy  (PI) 

Dr.  Wade  Davis  (Co-Investigator) 

Dr.  Tim  Hoffman  (Subcontract) 

Mr.  Matt  Gubin  (Graduate  student) 
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KEY  RESEARCH  ACCOMPLISHMENTS 


Outline  of  Statement  of  Work  (SOW)  (underlined  and  in  red  tasks  have  been  completed) 

Task  1 

Specific  A  im  1  :  Determine  HuR  targets  in  ER+  and  ER-  cell  lines  using  R  NAi  and 

immunoprecipitation. 

Experimental  Design: 

(1)  perform  transcriptomic  and  ribonomic  analyses  oil  ER+  and  ER-  cell  lines 

(2)  make  stable  ER  +  and  ER-  cell  lines  which  under  express  HuR  using  lentiviral  RNAi 

(3)  form  intersection  (Venn  diagram)  of  above  two  approaches  to  obtain  HuR  targets  in  ER+  and 

ER-  cell  lines 

(4)  validate  HuR  targets  using  Q-PCR 

Task  2 

Specific  Aim  2 :  Measure  growth  characteristics  of  ER+  and  ER-  cell  lines  w  ith  altered  HuR  levels 
ex  vivo  and  inject  into  nude  mice. 

Experimental  Design: 

(1)  make  stable  ER+  and  ER-  negative  cell  lines  which  over  express  epitope  tagged  HuR 

(2)  measure  growth  characteristics  of  altered  cell  lines  in  vitro 

(3)  inject  modified  ER+  and  ER-  cell  lines  into  athymic  nude  mice  and  monitor  tumor  growth 

(4)  remove  tumors  from  test  animals,  perform  transcriptomic  and  ribonomic  analyses  to  m  easure 

HuR  targets  using  Q-PCR,  Western  blotting  and  immunohistochemical  staining 

Task  3 

Specific  Aim  3:  Determine  whether  HuR  knockdown  by  RNAi  shrinks  established  tumors  in  nude 
mice. 

Experimental  Design: 

(1)  Establish  ER+  and  ER-  tumors  in  nude  mice 

(2)  Inject  lentivirus  expressing  HuR  RNAi,  measure  tumor  growth 

(3)  Remove  tumors  and  perform  transcriptomic  and  ribonomic  analyses  to  determine  effects  upon 
HuR  and  top  ten  targets  using  Q-PCR,  Western  and  immunohistochemical  staining 


•  Identification  of  distinct  sets  of  HuR  associated  mRNAs  from  ER-  and  ER+  breast  cancer  lines 

•  Development  of  stably  transfec  ted  ER+  and  ER-  breast  cancer  cell  lines  with  increas  ed  an  d 
decreased  HuR  levels 

•  Development  of  lentiviral  vectors  which  over-express  and  under-express  HuR 

•  Identification  and  validation  of  novel  HuR  associ  ated  mRNAs  from  ER-  and  ER+  breast  cancer 
(< cd9  and  calm2 ) 

•  Identification  of  novel  cancer  related  genes  which  are  HuR  targets 

•  HuR  diff  erentially  regu  lates  the  sa  me  mRNA  targe  ts  ( cd9  and  calm!)  depending  upon  cellular 
milieu 

•  Demonstration  that  changes  in  HuR  levels  in  ER-  breast  cancer  results  in  alterations  in  cell  cy  cle 
kinetics  and  altered  growth  rates 

•  Demonstration  that  HuR  over-expression  in  ER-  breast  cancer  increases  in  vitro  growth  rates  but 
in  orthotop  ic  m  odels  results  in  sig  nificant  decreases  (90  %)  in  turn  or  size  by  in  terfering  with 
angiogenesis 

•  Demonstration  that  HuR  over-expression  in  ER-  breast  cancer  resu  Its  in  s  ignificant  (75  %) 
reductions  in  angiogenesis 
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•  Knock-down  of  HuR  in  ER-  (MD  A-MB-23 1  cells)  using  lentiviral  shRNA  vectors  results  in 
significant  (95%)  reduction  in  tumor  size  in  orthotopic  models  of  breast  cancer 

Work  with  ER+  cell  lin  es  (MCF-7)  with  a  ltered  HuR  levels  was  ham  pered  by  their  inability  to 
grow  in  vitro.  Despite  several  attempts  at  re-m  aking  MCF-7  cells  with  higher  and  reduced  levels  of  HuR 
and  derving  stable  clones,  these  cells  did  not  appreciably  grow  in  vitro  and  therefore  we  were  unable  to 
transplant  into  athymic  nude  mice. 

We  made  lentiviral  vectors  which  over-expre  ssed  and  under-expressed  HuR  but  did  not  com  plete 
the  experim  ents  outlined  in  Task  3  (“  Determine  whether  HuR  k  nockdown  by  RNAi  shrinks 
established  tumors  in  nude  mice”),  due  to  tim  e  constraints.  Instead,  we  focused  and  followed  up  our 
surprising  findings  in  T  asks  1  and  2  that  HuR  over-e  xpression  reduced  turn  or  gr  owth  with  a  detailed 
investigation  into  mechanisms  of  tumor  reduction.  Our  putative  explan  ation  is  that  HuR  over-expression 
results  in  significant  reductions  in  angiogenesis  wh  ich  interf  eres  with  turn  or  growth.  The  mechanism 
involves  increases  in  T  SP- 1  (a  potent  angiogenic  i  nhibitor)  by  m  RNA  stabiliz  ation  and  a  concom  itant 
decrease  in  VEGF  (potent  pro-angiogenic  factor)  expression. 

As  per  instructions  requesting  complete  research  m  ethodologies,  we  have  in  eluded  two  m  anuscripts 
which  have  been  accepted  for  publication  based  upon  th  e  data  generated  above;  Ca  laluce  et  al  and  Gubin 
et  al.  these  two  manuscripts  are  included  as  attachments  in  the  Appendix.  These  manuscripts  describe  in 
detail  all  the  methodologies  and  research  results  in  detail. 


Milestones 

Year  1 

completed 

Year  2 

completed 

Year  3 

initiated 

Task  1 

completed 

Task  2 

mostly  completed 

Task  3 

initiated,  partially  completed,  not  finished 
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REPORTABLE  OUTCOMES 


(1)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  Advances  in  Breast  Cancer  Research:  Genetics,  Biology,  and  Clinical 
Applications  (San  Diego,  October  17-20,  2007) 

(2)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  AACR  Annual  Meeting  (San  Diego,  April  12-16,  2008) 

(3)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  Era  of  Hope  CDRMP  Meeting  (Baltimore,  June  25-28,  2008) 

(4)  Magee,  Joe,  Gubin,  Matt,  Atas  oy,  Ulus.  RNA  Immunopr  ecipitations  (RIP-on-Chip)  Identify  Unique 
Subsets  of  Genes  in  ER+  and  ER-  Breast  Cancer.  AACR  Annual  Meeting  (Denver,  April  18-22,  2009;  oral 
presentation) 

(5)  Calaluce,  Robert,  Gubin,  Matt,  Magee,  Joe,  Atasoy,  Ulus.  The  RN  A  binding  protein,  HuR,  regulates 
distinct  subsets  of  m  RNAs  in  estrogen  receptor  negative  (ER-)  and  estroge  n  rec  eptor  pos  itive  (ER+) 
breast  cancer.  AACR  Special  Symposium  on  Breast  Cancer  (San  Diego,  October  14-17,  2009;  poster 
presentation) 

(6)  Gubin,  Matt,  Calaluce,  Robert,  Magee,  Joe,  Casolaro,  Vincenzo,  Stellato,  Cristiana.  The  RNA  binding 
protein,  HuR,  coordinately  re  gulates  mRNAs  in  breast  cancer  and  C  D4+  Th2  polarization.  Keystone 
Symposium  (Keystone,  CO,  February  7-12,  2010;  poster  presentation). 

(7)  Gubin,  M,  Calaluce,  Bob,  Davis,  W  ade,  Magee,  Joe,  Strouse,  Connie,  Shaw,  Daniel,  Hoffm  an,  Tim, 
Rold,  Tammy,  Atasoy,  Ulus.  The  RNA  binding  protein  HuR  controls  a  ngiogenesis  in  triple  negative 
breast  cancer.  AACR  Annual  Meeting  (Washington,  DC,  April  15-21,  2010;  poster  presentation) 

(5)  Robert  Calaluce,  Matthew  M.  Gubin,  J.  Wade  Davis,  Joseph  D.  Magee,  Yuki  Kuwano,  Myriam 
Gorospe,  Ulus  Atasoy.  The  RNA  binding  protein,  HuR,  differentially  regulates  unique  subsets  of  mRNAs 
in  estrogen  receptor  negative  and  estrogen  receptor  positive  breast  cancer.  BMC  Cancer  2010,  10:126. 

(6)  Matthew  M.  Gubin,  Robert  Calaluce,  J.  W  ade  Davis,  Joseph  D.  Magee,  Connie  S. Strouse,  Daniel  P. 
Shaw,  Timothy  Hoffman,  Tammy  L.  Rold,  Ulus  Atasoy.  Over-expres  sion  of  the  RNA  binding  protein 
HuR  impairs  tumor  growth  in  triple  negative  breast  cancer  by  interferin  g  with  angiogenesis.  Cell  Cycle 
9:16,  1-11;  August  15,  2010. 
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CONCLUSIONS 

We  have  performed  RIP-Chip  (RNA  immunprecipitation  applied  to  microarrays,  (chips))  analysis 
from  ER-  and  ER+  breast  cancer  cell  lines  and  have  id  entitled  unique,  distinct  HuR  signature  subsets  of 
mRNAs  and  validated  novel  HuR  target  genes.  HuR  di  fferentially  regulates  the  same  mRNA  targets  in 
different  fashion  depend  ing  upon  cellular  m  ilieu.  Thes  e  m  ethods  m  ay  be  used  to  id  entify  novel  breast 
cancer  targets,  missed  by  traditional  gene  profiling  methods. 

Using  orthotopic  m  odels  of  breast  cancer  in  at  hymic  nude  m  ice,  we  have  m  ade  the  finding  that 
HuR  over  e  xpression  accelerates  MDA-MB-23 1  (ER-)  breast  cancer  growth  in  vitro,  but  surprisingly, 
results  in  a  significant  (90%)  reducti  on  in  turn  or  growth  as  com  pared  w  ith  controls  by  interf  ering  with 
angiogenesis.  The  putative  mechanisms  seem  to  be  an  anti-angiogenic  effect  in  that  HuR  over  expression 
increases  the  expression  of  throm  bospondin  1  (TSP-1)  a  potent  inhibitor  of  blood  vessel  form  ation  and 
simultaneously  decreases  VEGF.  In  parallel,  we  have  also  used  lentiviral  shRNA  to  knock-down  HuR  in 
ER-  breast  cancer;  when  inject  ed  into  athym  ic  nude  m  ice;  this  also  resulted  in  significant  reduction  in 
tumor  growth. 

“So  w  hat”:  These  m  ethods  will  enable  us  to  clarify  the  role  of  HuR  in  breas  t  cancer  as  well  as  th  e 
identification  of  novel  HuR  regulated  mRNA  targets  which  may  play  roles  in  breast  cancer.  Additionally, 
these  results  are  highly  significan  t  because  if  valid  ated,  HuR  m  ay  be  control  ling  a  posttran  scriptional 
operon  which  regulates  angiogenesis  in  ER-  tumor  fo  rmation.  Thus,  these  findings  m  ay  establish  a  ne  w 
paradigm  of  treatment  in  aggressive  triple  negative  breast  cancer  for  w  hich  currently  mono-therapies  do 
not  exist.  A  dditionally,  we  have  developed  rob  ust  lentiviral  vectors  which  can  over-express  and  under¬ 
express  HuR.  These  tools  can  be  used  to  infect  a  variety  of  cell  types. 
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APPENDIX 

(1)  See  attached  two  manuscripts,  Calaluce  et  al  (  BMC  Cancer  2010,  10:126)  and  Gubin  et  al.  (  Cell 
Cycle  9:16,  1-11;  August  15,  2010). 

http://www.ncbi.nlm.nih.gov/pubmed/203709 1 8 

http://www.ncbi.nlm.nih.gov/pubmed/20724828 

Please  see  end  of  Appendi  x  for  original  PDF  versions  of  these  tw  o  publicat  ions. 
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(2)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  Advances  in  Breast  Cancer  Research:  Genetics,  Biology,  and  Clinical 
Applications  (San  Diego,  October  17-20,  2007) 

Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by  the  RNA  Binding  Protein  HuR 
Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  University  of  Missouri,  Columbia,  Missouri. 

Due  to  the  poor  correlation  betw  een  steady  stat  e  m  RNA  levels  and  protein  products,  trad  itional 
microarray  analysis  may  miss  many  genes  which  are  regu  lated  primarily  at  the  level  of  mRNA  stability 
and  translation.  Posttranscriptional  gene  regulation  is  becoming  increasi  ngly  recognized  as  an  important 
form  of  cellular  control.  W  hereas  our  understanding  of  transc  riptional  regulation  is  advanced,  in  contrast 
this  area  remains  largely  unexplored.  Yet,  increa  singly,  the  importance  of  microRNAs  and  RNA  binding 
proteins  (RBPs)  is  now  beginning  to  be  better  appreciated.  We  study  the  elav  (embryonic  lethal  abnormal 
vision)  family  of  RBPs,  which  are  paraneoplastic  antigens,  over-express  ed  in  a  variety  of  m  alignancies, 
including  breast  cancer.  Antibodies  against  elav  family  members  are  believed  to  be  cancer-protective.  The 
elav  family  of  RBPs  binds  to  the  AU-rich  elements  (AREs)  found  in  the  3’  unt  ranslated  regions  (UTRs) 
of  many  early-response  genes,  including  proto-oncogenes  and  cell  cycle  regulators.  HuR,  the  ubiquitously 
expressed  family  member,  has  recently  been  described  to  play  a  role  in  cancer  progression.  HuR  stabilizes 
and  translationally  up  regulates  the  expression  of  its  target  mRNAs.  El  evated  levels  of  cytoplasm  ic  HuR 
directly  correlate  with  increased  invasiveness  of  m  alignancy  and  poor  prognos  is  for  m  any  cancers, 
including  those  of  the  breast.  HuR  has  been  descri  bed  to  control  the  expres  sion  of  genes  in  all  six 
different  areas  of  transform  ation  originally  descri  bed  by  Hanahan  and  W  einberg.  Hence,  it  has  been 
recently  suggested  that  HuR  may  serve  as  a  tumor  maintenance  gene  which  allows  for  cancers,  once  they 
are  established,  to  proliferate.  Therefore,  it  is  of  interest  to  discover  in  vivo  Hu  R  targets,  a  s  these  genes 
may  play  vital  roles  in  transform  ed  cells.  W  e  have  developed  methods  which  enable  us  to  identify  en 
masse,  in  vivo  targets  of  RBPs  such  as  HuR  fir  om  cell  lines  and  now  for  the  f  irst  time,  solid  tissues  as 
well.  W e  call  this  ribonom  ic  analysis,  in  which  we  i mmunoprecipitate  subsets  of  ribonuclear  particles 
(RNPs),  exact  the  m  RNA  and  then  apply  to  m  icroarrays.  T  his  technique,  in  which  no  cross  linking  is 
used,  allows  one  to  en  ma  sse  identify  different  in  vivo  mRNA  targets  which  different  RBPs  interact  with. 
Using  a  combination  of  HuR  i  mmunoprecipitation  and  RNA  interference,  we  have  been  able  to  identify 
HuR  in  viv  o  targ  ets  which  m  ay  play  im  portant  ro  les  in  tu  mor  m aintenance.  Furthermore,  we  pres  ent 
evidence  of  cross  talk  am  ongst  st  abilizer  and  de-stabilizer  RBPs  wh  ich  can  potentially  regulate  the 
expression  of  multiple  cancer  targets.  These  techniques  have  the  potential  to  better  delineate  gene  s  whose 
steady  state  mRNA  levels  may  not  significantly  change,  but  which  are  translationally  active  inside  cancer 
cells.  Potentially  the  identification  of  such  genes  m  ay  offer  novel  therapeutic  targets  to  inhibit  cancer 
growth. 


(3)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  AACR  Annual  Meeting  (San  Diego,  April  12-16,  2008) 

Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by  the  RNA  Binding  Protein  HuR 
Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  University  of  Missouri,  Columbia,  Missouri. 

Due  to  the  poor  correlation  betw  een  steady  stat  e  m  RNA  levels  and  protein  products,  trad  itional 
microarray  analysis  may  miss  many  genes  which  are  regu  lated  primarily  at  the  level  of  mRNA  stability 
and  translation.  Posttranscriptional  gene  regulation  is  becoming  increasi  ngly  recognized  as  an  important 
form  of  cellular  control.  W  hereas  our  understanding  of  transc  riptional  regulation  is  advanced,  in  contrast 
this  area  remains  largely  unexplored.  Yet,  increa  singly,  the  importance  of  microRNAs  and  RNA  binding 
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proteins  (RBPs)  is  now  beginning  to  be  better  appreciated.  We  study  the  elav  (embryonic  lethal  abnormal 
vision)  family  of  RBPs,  which  are  paraneoplastic  antigens,  over-express  ed  in  a  variety  of  m  alignancies, 
including  breast  cancer.  Antibodies  against  elav  family  members  are  believed  to  be  cancer-protective.  The 
elav  family  of  RBPs  binds  to  the  AU-rich  elements  (AREs)  found  in  the  3’  unt  ranslated  regions  (UTRs) 
of  many  early-response  genes,  including  proto-oncogenes  and  cell  cycle  regulators.  HuR,  the  ubiquitously 
expressed  family  member,  has  recently  been  described  to  play  a  role  in  cancer  progression.  HuR  stabilizes 
and  translationally  up  regulates  the  expression  of  its  target  mRNAs.  El  evated  levels  of  cytoplasm  ic  HuR 
directly  correlate  with  increased  invasiveness  of  m  alignancy  and  poor  prognos  is  for  m  any  cancers, 
including  those  of  the  breast.  HuR  has  been  descri  bed  to  control  the  expres  sion  of  genes  in  all  six 
different  areas  of  transform  ation  originally  descri  bed  by  Hanahan  and  W  einberg.  Hence,  it  has  been 
recently  suggested  that  HuR  may  serve  as  a  tumor  maintenance  gene  which  allows  for  cancers,  once  they 
are  established,  to  proliferate.  Therefore,  it  is  of  interest  to  discover  in  vivo  Hu  R  targets,  a  s  these  genes 
may  play  vital  roles  in  transform  ed  cells.  W  e  have  developed  methods  which  enable  us  to  identify  en 
masse,  in  vivo  targets  of  RBPs  such  as  HuR  fir  om  cell  lines  and  now  for  the  f  irst  time,  solid  tissues  as 
well.  W e  call  this  ribonom  ic  analysis,  in  which  we  i mmunoprecipitate  subsets  of  ribonuclear  particles 
(RNPs),  exact  the  m  RNA  and  then  apply  to  m  icroarrays.  T  his  technique,  in  which  no  cross  linking  is 
used,  allows  one  to  en  ma  sse  identify  different  in  vivo  mRNA  targets  which  different  RBPs  interact  with. 
Using  a  combination  of  HuR  i  mmunoprecipitation  and  RNA  interference,  we  have  been  able  to  identify 
HuR  in  viv  o  targ  ets  which  m  ay  play  im  portant  ro  les  in  tu  mor  m aintenance.  Furthermore,  we  pres  ent 
evidence  of  cross  talk  am  ongst  st  abilizer  and  de-stabilizer  RBPs  wh  ich  can  potentially  regulate  the 
expression  of  multiple  cancer  targets.  These  techniques  have  the  potential  to  better  delineate  gene  s  whose 
steady  state  mRNA  levels  may  not  significantly  change,  but  which  are  translationally  active  inside  cancer 
cells.  Potentially  the  identification  of  such  genes  m  ay  offer  novel  therapeutic  targets  to  inhibit  cancer 
growth. 
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(4)  Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by 
the  RNA  Binding  Protein  HuR.  Era  of  Hope  CDRMP  Meeting  (Baltimore,  June  25-28,  2008) 

Coordinate  Posttranscriptional  Regulation  of  Cancer  Genes  by  the  RNA  Binding  Protein  HuR 
Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus.  University  of  Missouri,  Columbia,  Missouri. 

Due  to  the  poor  correlation  betw  een  steady  stat  e  m  RNA  levels  and  protein  products,  trad  itional 
microarray  analysis  may  miss  many  genes  which  are  regu  lated  primarily  at  the  level  of  mRNA  stability 
and  translation.  Posttranscriptional  gene  regulation  is  becoming  increasi  ngly  recognized  as  an  important 
form  of  cellular  control.  W  hereas  our  understanding  of  transc  riptional  regulation  is  advanced,  in  contrast 
this  area  remains  largely  unexplored.  Yet,  increa  singly,  the  importance  of  microRNAs  and  RNA  binding 
proteins  (RBPs)  is  now  beginning  to  be  better  appreciated.  We  study  the  elav  (embryonic  lethal  abnormal 
vision)  family  of  RBPs,  which  are  paraneoplastic  antigens,  over-express  ed  in  a  variety  of  m  alignancies, 
including  breast  cancer.  Antibodies  against  elav  family  members  are  believed  to  be  cancer-protective.  The 
elav  family  of  RBPs  binds  to  the  AU-rich  elements  (AREs)  found  in  the  3’  unt  ranslated  regions  (UTRs) 
of  many  early-response  genes,  including  proto-oncogenes  and  cell  cycle  regulators.  HuR,  the  ubiquitously 
expressed  family  member,  has  recently  been  described  to  play  a  role  in  cancer  progression.  HuR  stabilizes 
and  translationally  up  regulates  the  expression  of  its  target  mRNAs.  El  evated  levels  of  cytoplasm  ic  HuR 
directly  correlate  with  increased  invasiveness  of  m  alignancy  and  poor  prognos  is  for  m  any  cancers, 
including  those  of  the  breast.  HuR  has  been  descri  bed  to  control  the  expres  sion  of  genes  in  all  six 
different  areas  of  transform  ation  originally  descri  bed  by  Hanahan  and  W  einberg.  Hence,  it  has  been 
recently  suggested  that  HuR  may  serve  as  a  tumor  maintenance  gene  which  allows  for  cancers,  once  they 
are  established,  to  proliferate.  Therefore,  it  is  of  interest  to  discover  in  vivo  Hu  R  targets,  as  these  genes 
may  play  vital  roles  in  transform  ed  cells.  W  e  have  developed  methods  which  enable  us  to  identify  en 
masse,  in  vivo  targets  of  RBPs  such  as  HuR  fir  om  cell  lines  and  now  for  the  f  irst  time,  solid  tissues  as 
well.  W e  call  this  ribonom  ic  analysis,  in  which  we  i mmunoprecipitate  subsets  of  ribonuclear  particles 
(RNPs),  exact  the  m  RNA  and  then  apply  to  m  icroarrays.  T  his  technique,  in  which  no  cross  linking  is 
used,  allows  one  to  en  ma  sse  identify  different  in  vivo  mRNA  targets  which  different  RBPs  interact  with. 
Using  a  combination  of  HuR  i  mmunoprecipitation  and  RNA  interference,  we  have  been  able  to  identify 
HuR  in  viv  o  targ  ets  which  m  ay  play  im  portant  ro  les  in  tu  mor  m aintenance.  Furthermore,  we  pres  ent 
evidence  of  cross  talk  am  ongst  st  abilizer  and  de-stabilizer  RBPs  wh  ich  can  potentially  regulate  the 
expression  of  multiple  cancer  targets.  These  techniques  have  the  potential  to  better  delineate  gene  s  whose 
steady  state  mRNA  levels  may  not  significantly  change,  but  which  are  translationally  active  inside  cancer 
cells.  Potentially  the  identification  of  such  genes  m  ay  offer  novel  therapeutic  targets  to  inhibit  cancer 
growth. 
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(5)  Magee,  Joe,  Gubin,  Matt,  Atas  oy,  Ulus.  RNA  Immunopr  ecipitations  (RIP-on-Chip)  Identify  Unique 
Subsets  of  Genes  in  ER+  and  ER-  Breast  Cancer.  AACR  Annual  Meeting  (Denver,  April  18-22,  2009;  oral 
presentation) 

RNA  Immunoprecipitations  (RIP-on-Chip)  Identify  Unique  Subsets  of  Genes  in  ER+  and  ER-  Breast 

Cancer 

Magee,  Joe,  Gubin,  Matt,  Atasoy,  Ulus,  University  of  Missouri,  Columbia,  Missouri. 

Due  to  the  poor  correlation  betw  een  steady  stat  e  m  RNA  levels  and  protein  products,  trad  itional 
microarray  analysis  may  miss  many  genes  which  are  regu  lated  primarily  at  the  level  of  mRNA  stability 
and  translation.  Posttranscriptional  gene  regulation  is  becoming  increasi  ngly  recognized  as  an  important 
form  of  cellula  r  contro  1.  The  im  portance  of  m  icroRNAs  and  RNA  binding  pr  oteins  (RBPs)  is  now 
beginning  to  be  better  appreciated.  We  study  the  elav  (embryonic  lethal  abnormal  vision)  family  of  RBPs, 
which  are  p  araneoplastic  antigens,  over-expressed  in  a  variety  of  m  alignancies,  including  breast  cancer. 
Antibodies  against  elav  fa  mily  members  are  believed  to  be  cancer-protective.  The  elav  family  of  RBPs 
binds  to  the  AU-rich  elem  ents  (AREs)  found  in  the  3  ’  untranslated  regions  (UTRs)  of  m  any  early- 
response  genes,  including  proto-oncogenes  and  cell  cy  cle  regulators.  HuR,  the  ubiquitously  expressed 
family  member,  has  been  described  to  play  a  role  in  cancer  progression.  HuR  stabilizes  and  translationally 
up  regulates  the  expression  of  its  target  m  RNAs.  Elevated  levels  of  cytoplasm  ic  HuR  directly  correlate 
with  increased  invasiveness  of  m  alignancy  and  poor  prognosis  for  m  any  cancers,  including  those  of  the 
breast.  Recently,  HuR  has  been  implicated  in  tarn  oxifen  resistance  in  breast  cancer.  HuR  has  been 
described  to  control  the  expression  of  genes  in  a  11  s  ix  dif  ferent  areas  of  tran  sfonnation  origin  ally 
described  by  Hanahan  and  W  einberg.  Hence,  it  has  been  suggested  that  HuR  may  serve  as  a  turn  or 
maintenance  gene  which  allows  for  cancers,  once  they  ar  e  established,  to  proliferat  e.  Therefore,  it  is  of 
interest  to  discover  in  vivo  HuR  targets,  as  these  genes  may  play  vital  roles  in  transformed  cells.  We  have 
developed  methods  which  enab  le  us  to  identify  en  m  asse,  in  vivo  ta  rgets  of  RBPs  such  as  HuR  from  cell 
lines  and  now  for  the  first  tim  e,  solid  tissues  as  well.  W  e  call  this  ribonomic  analysis  or  RNA 
immunoprecipitations  applied  to  microarrays,  RIP-on-Chip.  This  technique,  in  which  no  cross  linking  is 
used,  allows  one  to  en  ma  sse  identify  different  in  vivo  mRNA  targets  which  different  RBPs  interact  with. 
We  have  used  RIP-on-Chip  to  identify  distinct  subsets  of  HuR  associated  m  RNAs  in  MDA  MB23 1  and 
MCF-7  breast  cancer  cell  lines.  W  hereas  some  gene  targets  are  well  app  reciated  cancer  genes,  a  number 
do  not  have  clear  cut  associations  with  cancer  and  may  well  represent  novel  targets.  To  further  investigate 
the  role  of  HuR  in  triple  negative  breast  cancer,  we  have  over  expressed  HuR  in  MB23 1  cells,  which 
results  in  accelerated  growth.  Therefore,  we  cone  lude  that  HuR  pull  downs  can  identify  not  only  known 
but  also  novel  cancer  targets  and  RIP-on-Chip  analysis  identifies  distinct  subsets  of  cancer  relevant  genes. 
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Abstract 

Background:  The  discordance  between  steady-state  levels  of  mRNAs  and  protein  has  been  attributed  to 
posttranscriptional  control  mechanisms  affecting  mRNA  stability  and  translation.  Traditional  methods  of  genome 
wide  microarray  analysis,  profiling  steady-state  levels  of  mRNA,  may  miss  important  mRNA  targets  owing  to 
significant  posttranscriptional  gene  regulation  by  RNA  binding  proteins  (RBPs). 

Methods:  The  ribonomic  approach,  utilizing  RNA  immunoprecipitation  hybridized  to  microarray  (RIP-Chip), 
provides  global  identification  of  putative  endogenous  mRNA  targets  of  different  RBPs.  HuR  is  an  RBP  that  binds  to 
the  AU-rich  elements  (ARE)  of  labile  mRNAs,  such  as  proto-oncogenes,  facilitating  their  translation  into  protein.  HuR 
has  been  shown  to  play  a  role  in  cancer  progression  and  elevated  levels  of  cytoplasmic  HuR  directly  correlate  with 
increased  invasiveness  and  poor  prognosis  for  many  cancers,  including  those  of  the  breast.  HuR  has  been 
described  to  control  genes  in  several  of  the  acquired  capabilities  of  cancer  and  has  been  hypothesized  to  be  a 
tumor-maintenance  gene,  allowing  for  cancers  to  proliferate  once  they  are  established. 

Results:  We  used  HuR  RIP-Chip  as  a  comprehensive  and  systematic  method  to  survey  breast  cancer  target  genes 
in  both  MCF-7  (estrogen  receptor  positive,  ER+)  and  MDA-MB-231  (estrogen  receptor  negative,  ER-)  breast  cancer 
cell  lines.  We  identified  unique  subsets  of  HuR-associated  mRNAs  found  individually  or  in  both  cell  types.  Two 
novel  HuR  targets,  CD9  and  CALM2  mRNAs,  were  identified  and  validated  by  quantitative  RT-PCR  and  biotin  pull¬ 
down  analysis. 

Conclusion:  This  is  the  first  report  of  a  side-by-side  genome-wide  comparison  of  HuR-associated  targets  in  wild 
type  ER+  and  ER-  breast  cancer.  We  found  distinct,  differentially  expressed  subsets  of  cancer  related  genes  in  ER+ 
and  ER-  breast  cancer  cell  lines,  and  noted  that  the  differential  regulation  of  two  cancer-related  genes  by  HuR  was 
contingent  upon  the  cellular  environment. 


Background 

Over  the  past  decade  array  technologies  have  provided 
several  new  means  for  profiling  global  changes  in  gene 
expression.  The  power  of  DNA  microarrays  is  perhaps 
best  illustrated  in  the  way  it  has  been  used  to  differenti¬ 
ate  treatment  responses  in  patient  populations.  Indivi¬ 
dualized  and  targeted  therapy  for  several  tumors,  based 
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upon  underlying  differences  at  the  molecular  level 
among  gene  expression  profiles,  is  beginning  to  replace 
the  traditional  morphological-based  treatment  paradigm 
[1-3].  Genome  wide  microarray  analyses,  however,  are 
inherently  flawed  since  they  globally  profile  the  steady- 
state  levels  of  mRNA,  referred  to  as  the  transcriptome. 
Cellular  protein  expression  levels,  however,  do  not 
directly  correlate  with  steady-state  levels  of  mRNAs.  It 
is  well  accepted  in  the  RNA  field  that  there  is  a  poor 
correlation  between  steady-state  RNA  levels  and  protein. 
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This  discordance  has  been  attributed  to  posttranscrip- 
tional  control  mechanisms  affecting  mRNA  stability  and 
translation.  Steady-state  mRNA  levels  of  genes,  con¬ 
trolled  partially  or  totally  at  this  level,  may  be  mislead¬ 
ing.  Gygi  and  colleagues  have  shown  that  correlation 
between  mRNA  and  protein  levels  could  not  be  pre¬ 
dicted  from  only  mRNA  steady-state  levels  [4],  They 
observed  that  some  genes  had  the  same  mRNA  levels 
but  protein  levels  varied  more  than  20  fold.  Conversely, 
some  proteins  were  of  equal  expression  but  their  respec¬ 
tive  mRNA  level  varied  by  more  than  30-fold.  They  con¬ 
cluded  that  “transcript  levels  provide  little  predictive 
value  with  respect  to  the  extent  of  protein  expression” 
[4].  Additionally,  Idekar  and  colleagues  have  described 
similar  results  for  the  galactose  gene  [5]. 

Although  our  understanding  of  transcriptional  gene 
regulation  is  advanced,  posttranscriptional  gene  regula¬ 
tion  remains  largely  unexplored.  It  is  becoming  clear, 
however,  that  this  is  an  important  mode  of  gene  regula¬ 
tion,  especially  for  proinflammatory  genes.  These  genes 
appear  to  be  posttranscriptionally  regulated  by  RNA 
binding  proteins  (RBPs)  which  interact  with  AU-rich 
elements  (AREs)  in  the  3’  untranslated  region  (UTR)  of 
mRNAs.  Approximately  3,000  human  genes  contain 
AREs,  representing  8%  of  the  human  genome  [6].  Many 
of  these  genes  which  possess  AREs  are  in  areas  of  tran¬ 
sient  biological  responses,  including  cell  growth  and  dif¬ 
ferentiation,  immune  responses,  signal  transduction, 
transcriptional  and  translational  control,  hematopoiesis, 
apoptosis,  nutrient  transport,  and  metabolism  [6,7]. 

New  methodologies  have  provided  global  identifica¬ 
tion  of  in  vivo  mRNA  targets  of  different  RBPs.  One  of 
these,  termed  the  ribonomic  approach,  involves  the 
immunoprecipitation  of  ribonucleoprotein  complexes 
(RNPs)  with  antibodies  against  different  RBPs,  extrac¬ 
tion  of  mRNA,  and  hybridization  to  microarrays  [8-10]. 
This  approach,  also  referred  to  as  RIP-Chip,  enables 
investigators  to  identify  groups  of  posttranscriptionally 
regulated  mRNAs  coordinately  controlled  by  RBPs  dur¬ 
ing  various  biological  processes.  A  new  paradigm,  the 
posttranscriptional  operon  hypothesis,  has  been  devel¬ 
oped  which  states  that  RBPs  coordinately  regulate  the 
expression  of  biologically  related  molecules  [11,12].  This 
paradigm  is  being  confirmed  by  the  work  of  many  dif¬ 
ferent  laboratories  as  our  understanding  of  posttran¬ 
scriptional  regulation  broadens  and  putative  operons  are 
described  [8,13-17].  HuR  is  an  RBP  that  binds  to  AREs 
of  many  proto-oncogenes  and  labile  mRNAs.  It  has 
emerged  as  a  key  regulatory  factor  which  stabilizes  and 
translationally  enhances  its  targets  mRNAs,  and  affects 
their  transport  from  the  nucleus  to  the  cytoplasm 
[18-20].  HuR  belongs  to  the  ELAV  ( embryonic  lethal 
abnormal  vision)  family  found  in  mammalian  cells  con¬ 
taining  four  members:  HuR,  HuB,  HuC,  and  HuD.  HuR 


is  the  only  ubiquitously  expressed  member.  The  other 
family  members  are  found  primarily  in  the  central  ner¬ 
vous  system  and  gonadal  tissue  [18].  Many  HuR  targets 
are  cytokines,  chemokines,  and  other  early-response 
genes  [21,22]. 

Of  the  hallmarks  of  cancer  originally  described  by 
Hanahan  and  Weinberg,  HuR  has  been  demonstrated  to 
control  expression  of  genes  in  multiple  areas  of  malig¬ 
nant  transformation  [23],  Consequently,  HuR  has  been 
suggested  to  function  as  a  tumor  maintenance  gene,  per¬ 
missive  for  malignant  transformation,  tumor  growth,  and 
perhaps  metastasis  [24].  HuR  has  been  described  in  the 
literature  as  controlling  the  expression  of  many  cancer¬ 
relevant  genes,  including  those  that  encode  these  pro¬ 
teins:  Prothymosin-a,  Bcl-2,  Mcl-1,  SirTl,  TGF-|3,  MMP- 
9,  MTC-1,  uPA,  VEGF-a,  HIFl-a  and  cyclins  Al  (CCN 
Al),  B1  and  D1  [25-35].  Increased  levels  of  HuR  have 
been  associated  with  a  more  aggressive  breast  cancer  and 
a  worse  prognosis  [36-38].  Of  significance,  HuR  has  been 
described  as  posttranscriptionally  regulating  the  expres¬ 
sion  of  many  breast  cancer  relevant  genes  including 
those  that  encode  Glut-1,  ERa,  COX-2,  IL-8,  Cyclin  El, 
and  most  recently  BRCA-1  [36,39-44].  HuR  RIP-Chip 
analysis  has  recently  identified  Thrombospondin  1  as  a 
key  HuR  target  in  the  MCT-1  transformed  estrogen 
receptor  positive  (ER+)  cell  line  MCF-7  [45].  Its  interac¬ 
tions,  however,  are  complex  and,  at  times,  HuR  may 
interact  with  miRNAs  such  as  Let-7  to  translationally 
suppress  the  expression  of  C-MYC  mRNA  [46]. 

Since  HuR  has  been  described  as  regulating  the 
expression  of  many  cancer  relevant  genes,  we  asked 
whether  it  may  coordinately  regulate  breast  cancer 
genes  in  ER+  and  ER-  breast  cancer.  We  performed  a 
HuR  RIP-Chip  analysis  on  MDA-MB-231  (ER-)  and 
MCF-7  (ER+)  cell  lines  to  identify  cancer-relevant  genes, 
not  known  to  be  regulated  by  HuR,  and  potential  novel 
breast  cancer  targets.  Our  studies  indicated  that  HuR 
was  associated  with  unique  subsets  of  mRNAs  in  each 
cell  line  as  well  as  a  subset  of  HuR  associated  mRNA 
targets  common  to  both.  We  chose  two  cancer-asso¬ 
ciated  genes,  CD9  and  CALMODULIN  2  (CALM2), 
highly  expressed  in  both  cell  lines,  and  functionally  vali¬ 
dated  the  role  of  HuR  in  regulating  their  expression. 
Unexpectedly,  HuR  differentially  regulated  the  same  tar¬ 
get,  CD9,  in  both  cell  lines  in  an  opposite  manner. 
Moreover,  we  found  presumptive  differential  regulation 
of  CALM2  by  HuR,  as  HuR  interacted  only  with 
CALM2  mRNA,  but  not  with  family  members  CALM1 
and  CALM3  mRNAs.  We  discovered  that  HuR  interacts 
with  many  breast  cancer-relevant  genes  not  previously 
known  to  be  controlled  by  HuR,  and  target  genes  which 
have  not  been  shown  to  be  cancer  related.  This  latter 
category  may  indeed  represent  novel  cancer  genes  dis¬ 
covered  by  HuR  RIP-Chip  analysis. 
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Methods 

Cells  in  culture 

The  MDA-MB-231  (MB-231)  and  MCF-7  cell  lines  were 
obtained  from  American  Type  Culture  Collection  (Man¬ 
assas,  VA).  The  cell  lines  were  maintained  at  37°C  in  a 
humidified  atmosphere  of  95%  air  and  5%  C02.  MB-231 
cells  were  grown  in  RPMI  (GIBCO*,  Invitrogen™,  Carls¬ 
bad,  CA)  containing  10%  fetal  calf  serum  (Hyclone, 
Thermo  Fisher  Scientific,  Waltham,  MA),  0.5  mM  L- 
glutamine  (GIBCO8),  25  mg/ml  glucose  (Sigma-Aldrich), 
HEPES  (GIBCO*)  and  sodium  pyruvate  (GIBCO*). 
MCF-7  cells  were  grown  in  DMEM  (GIBCO*)  supple¬ 
mented  with  10%  fetal  calf  serum. 

HuR  Immunoprecipitations  (RIP-Chip) 

HuR  RIP-Chip  was  performed  as  previously  described 
[8,47,48].  Briefly,  lysates  were  prepared  from  exponen¬ 
tially  growing  MB-231  and  MCF-7  cells.  Equal  amounts 
of  protein  lysates  were  used  (100-300  pg).  HuR  mono¬ 
clonal  antibody  3A2  (made  in  our  laboratory  from  the 
3A2  hybridoma,  generously  provided  by  Dr.  loan  Steitz, 
Yale  University,  New  Haven,  CT),  or  isotype  control 
IgGl  (BD  Biosciences,  San  Jose,  CA),  were  pre-coated 
onto  Protein  A  Sepharose  beads  (PAS)  and  extensively 
washed.  Lysates  from  each  cell  initially  were  pre¬ 
absorbed  with  30  pg  of  IgGl  and  then  removed  by  addi¬ 
tion  of  PAS  beads.  Individual  pull-downs  were  per¬ 
formed  at  4°C  for  only  1-2  hr  to  minimize  potential  re¬ 
assortment  of  mRNAs. 

RNA  amplification 

The  entire  amount  of  recovered  RNA  per  immunopreci- 
pitation  was  amplified  using  the  WT-Ovation™  Pico 
RNA  Amplification  System  protocol  (NuGen,  San  Car¬ 
los,  CA).  Forty  ng  of  total  RNA  was  used  as  starting 
material  to  generate  at  least  6  pg  of  cDNA.  Amplified 
cDNA  was  purified  using  Zymo  Research  Clean  and 
ConcentratorT"-25  (Zymo  Research,  Orange,  CA).  Three 
pg  of  amplified  and  purified  cDNA  was  incubated  at  50° 
C  for  30  minutes  with  5  pi  of  UNG  buffer  and  5  pi 
UNG  enzyme  and  60  minutes  with  5  pi  labeling  buffer 
and  5  pi  ARP  (biotin)  solution  as  described  in  NuGen’s 
labeling  protocol  for  the  Illumina  BeadArray  platform. 
All  samples  (total  RNA,  amplified  cDNA,  and  biotin 
labeled  amplified  cDNA)  were  quantitated  using  a 
Nanodrop™  (Thermo  Fisher  Scientific,  Waltham,  MA) 
spectrophotometer.  RNA  quality  and  integrity  were 
assessed  on  selected  samples  with  the  Experion™  auto¬ 
mated  electrophoresis  system  (Bio-Rad,  Hercules,  CA). 

Microarray 

Biotin-labeled,  amplified  cDNA  (1.5  pg)  was  hybridized 
to  a  Sentrix*  Human-6  v.2  Whole  Genome  Expression 


BeadChips  (Sentrix  Human  WG-6;  Illumina,  San  Diego, 
CA).  Each  chip  tested  6  samples  and  contained  47,293 
gene  targets,  representing  18,025  distinct  RefSeq  genes 
that  are  not  pseudogenes.  A  total  of  3  chips  were  used 
for  this  experiment.  The  chips  were  hybridized  at  48°C 
for  20  hr  in  the  hybridization  buffer  provided  by  the 
manufacturer.  After  hybridization,  the  chips  were 
washed  and  stained  with  streptavidin-C3.  The  chips 
were  scanned  on  the  BeadArray  Reader,  as  described  by 
Illumina  at  http://www.illumina.com.  The  Illumina 
BeadStudio  software  was  used  to  assess  fluorescent 
hybridization  signals. 

Quantitative  RT-PCR 

Selected  genes  were  validated  by  quantitative  RT-PCR. 
Briefly,  cDNA  was  generated  from  the  same  samples  as 
previously  described  for  the  microarray  experiments 
using  10  ng  total  RNA  and  the  Superscript™  III  Plati¬ 
num*  Two-Step  qRT-PCR  Kit  with  SYBR*  Green  (Invi¬ 
trogen  Carlsbad,  CA).  RT-PCR  was  performed  on  the 
StepOne™  Real-Time  PCR  System  (Applied  Biosystems, 
Foster  City,  CA).  Each  sample  was  run  in  triplicate  for 
these  genes  and  the  cDNA  was  divided  equally  per  reac¬ 
tion  in  a  20  pi  volume.  The  PCR  conditions  were:  50°C 
for  2  minutes  and  95°C  for  2  minutes,  followed  by  40 
cycles  of  95°C  for  15  seconds  alternating  with  60°C  for 
30  seconds.  Melting  curve  analysis  was  performed  on 
every  reaction  to  confirm  a  single  amplicon.  For  each 
cell  line,  differences  in  gene  expression  were  determined 
using  the  equation  2'AACt,  where  the  Ct  value  for  either 
the  HuR  or  IgG  IP  was  subtracted  from  the  Ct  value  of 
the  GAPDH  control  to  yield  the  A Ct  value.  For  each  cell 
line,  the  A Ct  value  for  the  HuR  and  IgG  IP  were  com¬ 
puted  in  triplicate  and  averaged  to  give  one  A  A  Ct  value 
per  sample.  Primers  used: 

Human  RT  GAPDH  Forward  5’  AGCCTCAAGAT- 
CATCAGCAATGCC  3’ 

Reverse  5’  TGTGGTCATGAGTCCTTCCACGAT  3’ 
Human  RT  HuR  Forward  5’  ATGAAGACCA- 
CATGGCCG AAG ACT  3’ 

Reverse  5’  AG  T  T  C  AC  A  A  AG  CC  AT  AG  CCC  A  AG  C  3’ 
Human  RT  CD9  Forward  5’  TCAGACCAAGAG- 
CAT CTTCG AGCA  3’ 

Reverse  5’  ACCAAGAGGAAGCCGAAGAACAGT  3’ 
Human  RT  CALM2  Forward  5’  CT G ACCAACT - 
G ACT G AAG AGCAG A  3’ 

Reverse  5’  TTCTGTGGGATTCTGCCCAAGAG  3’ 

Cloning  strategy  of  HA  HuR 

Hemagglutinin  (HA)  tagged  human  HuR  was  cloned 
into  the  Nhel  and  Xhol  sites  of  the  pZeoSV2  (-)  vector 
(Invitrogen).  The  plasmids  were  sequenced  in  both 
directions  to  confirm  identity.  Cells  were  transfected 
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with  either  pZeo  HA  HuR  or  pZeo  empty  vector  using 
Lipofectamine  2000  (Invitrogen).  After  five  days  trans¬ 
fected  media  was  removed  and  replaced  with  fresh  med¬ 
ium  containing  200  pg/ml  of  Zeocin  antibiotic 
(Invitrogen).  Cells  were  selected  for  a  ten  day  period. 
After  ten  days,  the  selected  cells  were  maintained  in  50 
pg/ml  of  Zeocin  to  maintain  pZeo  HA  HuR  and  empty 
vector  expression.  No  viable  cells  remained  in  the 
untransfected  well.  Cells  were  then  cloned  by  limiting 
dilution. 

Lentiviral  RNAi  HuR  knock-down 

In  order  to  knockdown  HuR,  PSICOOLIGOMAKER 
vl.5  http://web.mit.edu/ccr/labs/jacks/was  used  to 
identify  optimal  shRNAs  sequences  to  HuR.  We  tested 
multiple  sequences  and  chose  GGATCCTCTGGCA- 
GATGT,  identified  and  designated  shRNA  H760. 
Annealed  sense  and  antisense  DNA  (Integrated  DNA 
Technologies,  Inc,  IDT,  Coralville,  IA),  along  with  stem 
loops  to  create  hairpin,  were  cloned  into  the  Hpal  and 
Xhol  restriction  sites  in  the  Lentilox  pll3.7  vector 
(ATCC).  After  sequence  verification,  lentivirus  was 
packaged  in  293FT  cells  using  ViraPower™  Lentiviral 
Expression  Systems  (Invitrogen)  following  manufac¬ 
turer’s  protocol.  Both  MB-231  and  MCF-7  cells  were 
seeded  at  a  density  of  100,000  cells  in  100  mm  tissue 
culture  plates  with  10  ml  of  media.  The  following  day 
lentivirus,  expressing  either  GFP  and  no  shRNA  (empty 
lentilox  control)  or  GFP  and  HuR  shRNA  H760,  was 
added  at  a  multiplicity  of  infection  (MOI)  of  10  along 
with  polybrene  (8  pg/ml)  (Sigma-Aldrich  Corp,  St. 
Louis,  MO).  After  five  days,  cells  were  harvested  by 
trypsinization  and  sorted  for  GFP  expression  using  BD 
FACSDiva  (BD  Bioscience).  Cells  were  cloned  by  limit¬ 
ing  dilution  and  GFP  expression  was  assessed  using 
FACScan  (BD  Bioscience)  and  CellQuest  software  (BD 
Bioscience).  GFP  expression  was  >98%  and  indicated 
homogenous  cell  population. 

SDS-PAGE  and  Western  Blot  Analysis 

Western  analysis  was  performed  as  described  previously 
with  slight  modifications  [47].  Briefly,  cells  were  har¬ 
vested  and  lysed  in  triple-detergent  RIPA  buffer  with 
protease  inhibitor  cocktail  (Roche,  Pleasanton,  CA).  For 
nuclear  and  cytoplasmic  fractionation,  the  NE-PER  kit 
was  used  (Pierce,  Rockford,  IL).  Protein  quantity  was 
determined  by  Bradford  Assay.  Forty  pg  of  protein  was 
electrophoresed  on  a  12%  SDS-polyacrylamide  gel  and 
transferred  to  a  nitrocellulose  membrane.  The  mem¬ 
brane  was  blocked  with  5%  nonfat  milk  at  room  tem¬ 
perature  for  1  hr  and  incubated  with  anti-P -tubulin  (1 
pg/ml,  Sigma-Aldrich)  at  4°C  overnight.  After  washing, 
the  membrane  was  incubated  with  monoclonal  anti- 
HuR  clone  3A2  antibody  (1  pg/ml)  at  room  temperature 


for  1  hr,  or  anti-CD9  antibody  (1:100)  (Santa  Cruz  Bio¬ 
technology,  Inc.,  Santa  Cruz,  CA)  at  4°C  overnight.  The 
secondary  antibody  used  was  sheep  anti-mouse  Ig  horse 
radish  peroxidase  (1:4000)  (GE  Healthcare,  Piscataway, 
NJ),  incubated  at  room  temperature  for  1  hr.  Specific 
proteins  were  detected  using  chemiluminescence  (GE 
Healthcare).  HuR  knock-down  was  determined  to  be 
>90%  using  Bio-Rad’s  Quantity  One  software  (Bio-Rad) 
normalizing  to  P -tubulin,  and  HuR  over-expression  was 
quantitated  in  a  similar  manner. 

Biotin  Pull-down 

Biotinylated  transcripts  were  synthesized  using  cDNA 
that  was  prepared  from  MB-231  cells.  Templates  were 
prepared  using  forward  primers  that  contained  the  T7 
RNA  polymerase  promoter  sequence  (CCAAGCTTC- 
TAATACGACTCACTATAGGGAGA  [T7]).  Primers 
used  for  the  preparation  of  biotinylated  transcripts  span¬ 
ning  the  CD9  CR,  and  3’UTR  (NM_001769)  and 
CALM2  CR  and  3’UTR  (NM_001743.3)  were  as  follows: 

CD9  CR  118-560:  [T7]  TCAAAGGAGGCAC- 
CAAGTGCAT and  AACGCATAGTGGATGGCTTTCA 

CD93’UTR798-1231:  [T7]  AGTCAGCTTACATC- 
CCTGAGCA  and  GACATTGTCATAATTTTTTAT- 
TATGTATC 

CALM2  CR  72-515:  [T7]  GCTGACCAACTGACT- 
GAAGA  and  CTTTGCTGTCATCATTTGTACAAA 

CALM2  3’UTR  518-1128:  [T7]  AGACCTTGTACA- 
GAATGTGTTAA  and  GGGTAAATTGTAATTTTT- 
TTATTGGAA 

GAPDH  3’UTR:  [T7]  CCT CAACGACCACTTTGTCA 
and  GGTT G AGCACAGGG  TACTTTATT 

The  PCR-amplified  fragments  were  purified  and  used 
as  templates  for  in  vitro  synthesis  of  the  corresponding 
biotinylated  RNAs  by  MAXIscript  kit  (Ambion®,  Applied 
Biosystems).  Biotin  pull-down  assays  were  performed  by 
incubating  40  pg  of  MB-231  cell  lysates  with  equimolar 
of  biotinylated  transcripts  for  1  hr  at  room  temperature. 
The  complexes  were  isolated  using  paramagnetic  strep- 
tavidin-conjugated  Dynabeads  (Dynal®,  Invitrogen),  and 
bound  proteins  in  the  pull-down  material  were  analyzed 
by  Western  blotting  using  an  antibody  recognizing  HuR 
(Santa  Cruz).  After  secondary- antibody  incubations,  the 
signals  were  visualized  by  chemiluminescence  (Amer- 
sham  Biosciences,  GE  Healthcare). 

Statistical  Analysis  of  Microarray  Data 

Analysis  of  microarray  gene  expression  data  was  primar¬ 
ily  performed  using  the  Linear  Models  for  Microarray 
Data  (limma)  package  [49]  and  the  lumi  package  [50], 
available  through  the  Bioconductor  project  [51]  for  use 
with  R  statistical  software  [52].  After  data  pre-processing 
was  completed  (Appendix),  the  statistical  analysis  was 
performed  using  moderated  t-statistics  applied  to  the 
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log-transformed  (base  2)  normalized  intensity  for  each 
gene  using  an  Empirical  Bayes  approach  [53].  Three 
contrasts  of  interest  were  computed  and  tested.  The 
first  was  the  difference  between  HuR  pull-down  and  IgG 
background  for  the  MB-231  cell  line.  Genes  which 
exhibited  significantly  greater  expression  in  the  pull¬ 
down  were  considered  to  be  in  the  HuR  pellet  for  the 
MB-231  cell  line.  The  second  contrast  was  similar  to 
the  first,  but  for  the  MCF-7  cell  line.  The  third  and 
most  important  contrast  was  the  difference  between  the 
first  and  second  contrast,  and  can  be  viewed  as  a  test  of 
statistical  interaction  between  HuR  and  cell  line.  For  a 
given  gene,  this  term  can  be  interpreted  as  reflection  of 
the  synergistic  relationship  between  HuR  and  estrogen 
in  breast  cancer.  Adjustment  for  multiple  testing  was 
made  using  the  false  discovery  rate  (FDR)  method  of 
Benjamini  and  Hochberg  [54]  with  an  FDR  of  10%  as 
our  cutoff  for  declaring  significance.  To  facilitate  inter¬ 
pretation,  log  fold  changes  were  transformed  back  to 
fold  change  on  the  data  scale  (fluorescent  intensity). 

Gene  ontology  (GO)  analyses  were  carried  out  on  the 
list  of  significant  genes  based  on  the  third  contrast 
described  above.  The  purpose  of  the  analyses  was  to  test 
the  association  between  Gene  Ontology  Consortium 
categories  [55]  and  differentially  expressed  HuR  pellet 
genes  between  MB-231  and  MCF-7.  Using  our  defined 
gene  universe  (Appendix),  GOstats  [56]  was  used  to 
carry  out  conditional  hypergeometric  tests.  These  tests 
exploit  the  hierarchical  nature  of  the  relationships 
among  the  GO  terms  for  conditioning  [57].  We  carried 
out  GO  analyses  for  over-representation  of  biological 
process  (BP),  molecular  function  (MF),  and  cellular 
component  (CC)  ontologies,  and  computed  the  nominal 
hypergeometric  probability  for  each  GO  category.  These 
results  were  used  to  assess  whether  the  number  of 
selected  genes  associated  with  a  given  term  was  larger 
than  expected  under  the  null  hypothesis,  and  a  /7-value 
cutoff  of  0.01  was  used.  GO  categories  containing  less 
than  10  genes  from  our  gene  universe  were  not  consid¬ 
ered  to  be  reliable  indicators,  and  are  not  reported. 

Results 

HuR  immunoprecipitation  from  ER+  and  ER-  breast  cancer 
cell  lines 

We  first  determined  HuR  protein  expression  levels  in 
breast  cancer  cell  lines.  HuR  is  expressed  in  both  the 
ER-  and  the  ER+  cell  lines,  MB-231  and  MCF-7,  respec¬ 
tively  (Figure  1A).  RNA  immunoprecipitation,  using 
HuR  monoclonal  antibody  3A2,  recovered  HuR  (Figure 
1A)  and  revealed,  by  quantitative  RT-PCR,  a  significant 
enrichment  of  up  to  fifteen  fold  for  a  known  HuR  tar¬ 
get,  /3-ACTIN  mRNA,  as  compared  to  isotype  control 
(IgGl)  and  normalized  to  a  non-target,  GAPDH  mRNA 
(Figure  IB).  These  data  showed  that  HuR  RIP 


specifically  immunoprecipitate  HuR  protein  and  asso¬ 
ciated  mRNAs,  though  absolute  quantitative  conclusions 
cannot  be  drawn  since  different  amounts  of  lysates  were 
used  and  efficiency  of  immunoprecipitation  from  differ¬ 
ent  cell  lines  may  differ. 

RIP-Chip  from  ER+  and  ER-  breast  cancer  cell  lines 
identifies  unique  sets  of  associated  mRNAs 

RIP-Chip  was  performed  on  cytoplasmic  lysates  from 
both  breast  cancer  cell  lines  with  HuR  antibody  and  iso¬ 
type  control  in  order  to  determine  HuR  associated 
mRNAs.  Each  immunoprecipitation  was  done  at  least 
three  independent  times  with  matching  controls.  Signals 
from  isotype  control  were  subtracted  out.  Recovered 
mRNA  was  amplified  and  hybridized  to  Illumina  Sentrix 
Human  arrays  consisting  of  47,000  genes.  Figure  2 
represents  a  composite  array  generated  by  combining 
hybridizations  to  twelve  different  arrays  (log2  scale). 
Three  groups  of  HuR-associated  target  genes  were  iden¬ 
tified:  MB-231  targets  in  the  left  upper  quadrant;  both 
MB-231  and  MCF-7  targets  in  the  right  upper  quadrant; 
MCF-7  targets  in  the  right  lower  quadrant.  As  expected, 
most  of  the  mRNAs  did  not  associate  with  HuR  and 
were  located  in  the  lower  left  quadrant.  There  were  395 
and  64  annotated  genes,  at  least  2  fold  or  more 
enriched,  associated  with  either  MB-231  or  MCF-7  cells, 
respectively,  and  182  genes  associated  with  both  cell 
lines.  A  complete  list  can  be  found  in  Additional  File  1, 
Figure  S2.  The  raw  data  files  are  available  in  the  NCBI 
database  at  the  following  link;  http://www.ncbi.nlm.nih. 
gov/geo/query/acc.cgi?token=pdsnrqmiawukql- 
m&acc=GSE17820,  NCBI  Accession  number 
GSE17820).  These  genes  generally  fell  into  three  groups. 
Group  1  consisted  of  cancer-associated  genes  which 
were  known  HuR  targets,  such  as  PTMA  mRNA.  Group 
2  consisted  of  genes  which  played  a  role  in  cancer  but 
were  not  known  to  be  HuR  targets.  Group  3  consisted 
of  genes  with  an  unknown  function  in  cancer,  but 
which  may  be  regulated  by  HuR.  These  data  revealed 
that  HuR  was  associated  with  distinct  subsets  of 
mRNAs  in  ER+  and  ER-  breast  cancer  cells. 

Gene  Ontology  (GO)  analyses  of  differentially  expressed 
significant  genes  between  ER+  and  ER-  cells  were  cate¬ 
gorized  into  Biological  Process  (BP),  Cellular  Compo¬ 
nent  (CC),  and  Molecular  Function  (MF).  GO  analyses 
allow  for  the  identification  of  gene  families  that  may 
play  significant  roles  related  to  these  categories  in 
expression  profiles.  Most  of  the  differentially  expressed 
genes  (155)  were  found  to  be  more  abundant  than 
expected  in  14  BP  categories  (Figure  3A).  Three  MF 
categories  consisted  of  100  genes  with  most  of  these 
(83)  related  to  protein  binding  and  transcription  activa¬ 
tor  activity.  The  CC  categories  contained  the  least  (34) 
and  were  primarily  associated  with  the  Golgi  apparatus. 
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Figure  1  Immunoprecipitation  and  RIP  in  MB-231  and  MCF-7  breast  cancer  cells  Immunoprecipitations  were  performed  from  MB-231  or 
MCF-7  cell  lysates  using  anti-HuR  monoclonal  antibody  (3A2)  and  IgGI  isotype  control.  A.  IP  Western  of  HuR  revealed  expected  size  band  as 
detected  by  3A2.  Panel  on  right  reveals  amounts  of  HuR  in  lysates  used  from  both  cell  lines.  B.  Verification  by  quantitative  RT-PCR  showed 
fifteen  and  eleven  fold  enrichments  of  B-ACTIN,  a  known  HuR  target,  in  the  3A2  IPs  from  MB231  and  MCF-7,  respectively.  All  AACT  values  were 
normalized  to  GAPDH.  Experiments  were  done  in  duplicate  (n  =  2). 


o 

</> 

3 

C 


co 


co 

cvi 

CD 


MB-231 


MB-231  &  MCF-7 


MCF-7 


MCF-7  3A2  minus  IgG 


Figure  2  HuR  RIP-CHIP  identifies  distinct  genetic  profiles  in  ER+  and  ER-  breast  cancer  cells  HuR  immunoprecipitations  were  performed 
from  MB-231  or  MCF-7  cell  lysates  using  HuR  antibody  and  IgGI  isotype  control  hybridized  to  lllumina  Sentrix  arrays  (47,000  genes).  Control 
signals  were  subtracted.  Results  represent  cumulative  data  from  12  different  arrays.  Experiments  were  done  in  triplicate  (n  =  3)  for  each  cell  line 
with  matching  controls.  Scales  are  log2. 


Calaluce  et  al.  BMC  Cancer  2010,  10:126 
http://www.biomedcentral.com/1471-2407/10/126 


Page  7  of  15 


Epithelial  ceil  differentiation 
Homi  one  metabolisn 
Regulation  of  biological  processes 
Blood  vessel  morphogenesis 
Anatomical  Sructure  formation 
Vasculature  development 
Nucleie  add  metabolism 
Macromolecule  biosynthesis 
Regulation  of  gene  expression 
Regulation  of  metabolic  pocestes 
Transcription  regulation.  DNA  dependent 
Transcription  regulation  from  RNA  poi  II  promoter 
Regulation  of  Cellular  Processes 
Signal  Transduction 


7  9 

Odds  Ratio 


B 


/ 


EOF  GAT  A3 

GMCSt  CAV1 

FOS  EREG 

Vf  VC  FOAAt 

CYGUHA 

CVCUHBt 

crcuHDt 


EREG 

CAV1 

NCOA3 


EREG 

NCOA3 

AFF4 


Survival  of 
Apoptotlc  Damage 


\ 


X 


TGF-d 

GALECm 


IL8 

RHOB 

CAV1 

RUNX1 


Elevated 

Local 

Angiogenesis 


Evasion  of 
Immune 
Recognition 


Figure  3  GO  Classification  of  genes  found  by  RIP  CHIP  of  potential  HuR  targets  and  their  relationship  to  the  Acquired  Capabilities  of 
Cancer  Model.  A.  Differentially  expressed  genes  which  are  more  represented  in  the  Biological  Processes  (BP)  GO  category  than  expected.  B. 
Original  representation  showing  subsets  of  transcripts  found  to  be  targets  of  association  with  HuR  (normal  type).  New  transcripts  found  in  this 
study  with  RIP-Chip  (bold  type).  Enhanced  expression  upon  binding  to  HuR  influences  several  of  the  acquired  capabilities  of  cancer  cells 
described  by  Hanahan  and  Weinberg  [23,24], 


For  the  complete  GO  analyses  see  Additional  File  2, 
Table  SI.  In  Table  SI  we  list  the  top  FluR  associated 
mRNAs  in  the  different  categories  which  were  approxi¬ 
mately  5  fold  enriched  or  greater.  As  can  be  seen  in  Fig¬ 
ure  3B,  a  partial  listing  of  some  of  these  genes  (in  bold) 
are  candidate  members  to  multiple  areas  of  cancer  con¬ 
trol,  as  suggested  by  Hanahan  and  Weinberg.  We  note 
that  though  P -ACTIN  mRNA  was  amongst  the  most 
abundant  of  HuR-associated  mRNAs  in  MCF-7  cells,  P- 
ACTIN  mRNA  levels  were  only  3.93-fold  higher  in  HuR 
IP  compared  to  IgG  IP  in  MB-231  cells.  Therefore, 
since  this  was  less  than  the  5-fold  cut-off  we  employed 
for  Table  SI,  it  is  not  listed.  Thus,  these  results  may 
have  identified  novel  HuR-controlled  genes  which  may 


play  roles  in  breast  carcinogenesis  in  a  cancer  subtype- 
specific  fashion. 

Validation  of  HuR  targets  CD9  and  CALM2  by  real-time 
PCR  and  biotin  pull-down  analyses 

In  order  to  validate  HuR  binding  to  genes  identified  in  Fig¬ 
ure  2,  we  chose  two  known  cancer  associated  genes,  CD9 
and  CALM2,  highly  expressed  in  both  cell  lines.  Two  inde¬ 
pendent  approaches  confirmed  the  physical  interaction 
between  HuR,  CD9  and  CALM2  mRNAs.  Precipitated 
mRNA  from  the  RIP-Chip  experiments  were  analyzed  by 
quantitative  RT-PCR.  Both  CD9  and  CALM2  mRNAs 
were  enriched  in  the  HuR  RIP  by  as  much  as  160  fold 
(Figures  4A  and  4B),  but  not  the  isotype  control  IP.  We 
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Figure  4  Validation  of  target  CALM2  and  CD9  mRNAs  by  quantitative  RT-PCR.  Quantitative  RT-PCR  using  cell  lysates,  HuR  antibody,  and 
IgGI  from  RIP-CHIP  analysis  confirmed  results  identifying  CALM2  mRNA  (A)  and  CD9  mRNA  (B)  as  HuR  targets.  Change  in  gene  expression  is 
represented  as  fold  increase  in  HuR  immunoprecipitation  as  compared  to  IgGI.  GAPDH  mRNA  was  used  as  an  endogenous  control.  Error  bars 
represent  SEM.  p  value  is  <  0.005.  Experiments  were  done  in  triplicate  (n  =  3). 


further  confirmed  HuR  binding  to  CD9  and  CALM2 
mRNAs  by  biotin  pull-downs.  The  relevant  portion  of  the 
mRNA  was  transcribed  with  biotin  tags  and  incubated 
with  lysates  from  the  two  cell  lines  to  probe  for  interac¬ 
tions  with  protein.  The  mixture  was  then  separated  by 
pull-down  using  streptavidin-coated  beads  and  HuR  levels 
were  analyzed  by  Western  blot  analysis.  As  seen  in  Figure 
5,  HuR  specifically  interacts  with  CD9  and  CALM2 
mRNAs  in  the  3’UTR  regions,  but  not  within  the  coding 
region  (CR)  or  with  a  control  biotinylated  RNA  corre¬ 
sponding  to  the  3’UTR  of  the  housekeeping  control 
GAPDH  mRNA,  which  is  not  a  target  of  HuR. 


HuR  differentially  regulates  CD9  and  CALM2  in  MB-231 
and  MCF-7  cell  lines 

To  gain  insight  into  the  biological  effects  of  these  asso¬ 
ciations,  we  studied  the  consequences  of  stably  increas¬ 
ing  or  decreasing  HuR  abundance.  Individual  MB-231 
clones  which  over-  and  under-express  HuR  were  estab¬ 
lished  by  limiting  dilution  (Figure  6A  and  6B).  MB-231 
cells  over-expressed  HuR  by  about  140%  (Figure  6A). 
HuR  knock-down  using  lentiviral  shRNA  resulted  in 
~95%  reduction  in  HuR  expression  (Figure  6B).  Surpris¬ 
ingly,  over-expression  of  HuR  in  MB-231  cells  caused 
decreases  in  both  CD9  protein  and  mRNA  levels 
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Figure  5  Biotin  Pull-down  of  CD9  and  CALM2.  A.  Scheme  of  Coding  region  (CR)  and  3'UTR  fragments  for  biotin  pull-down  assay.  The 
sequences  were  obtained  from  Entrez  data  base.  CR  and  3'UTR  fragments  selected  for  amplification  by  PCR  are  as  noted.  B.  1%  agarose  gel 
electrophoresis  showing  PCR  amplified  products  of  the  coding  regions  and  3'UTR's  for  CD9  (442  bp  and  432  bp,  respectively)  and  CALM2  (443 
bp  and  610  bp,  respectively).  C.  Biotin  pull-down  assay  using  lysates  prepared  from  MB-231  cells.  The  binding  of  HuR  to  biotinylated  3'UTR 
transcripts  from  both  CD9  and  CALM2  mRNAs  was  specific.  HuR  did  not  bind  a  biotinylated  control  (GAPDH  3'UTR);  and  did  not  bind  to 
biotinylated  transcripts  spanning  the  CR  of  CD9  or  CALM2.  Experiments  were  done  in  duplicate  (n  =  2). 


(Figures  6C  and  6D).  HuR  knock-down,  however, 
resulted  in  increases  in  both  CD9  mRNA  and  protein 
levels  (Figures  6C  and  6E).  This  is  the  opposite  of  what 
we  predicted,  since  HuR  is  generally  regarded  as  a  stabi¬ 
lizer  of  mRNA.  In  contrast,  over-expression  of  HuR  in 
MB-231  cells  did  not  significantly  alter  the  levels  of 
CALM2  mRNA  (Figure  6D).  Figure  6F  depicts  a  graphi¬ 
cal  analysis  which  reveals  that  HuR  over-expression 
decreases  both  CD9  mRNA  and  protein  levels,  as  com¬ 
pared  to  controls  (dashed  line  set  at  100%).  Whereas, 
HuR  shRNA  knock-down  results  in  increases  in  both 
CD9  mRNA  and  protein  levels  above  control  levels. 

We  performed  similar  analyses  with  MCF-7  cells, 
though  the  over-expression  levels  of  HA  HuR  were  only 
approximately  10%,  since  this  was  a  pooled  population 
(we  have  been  unable  to  obtain  MCF-7  clones  which 


over-express  HuR).  In  contrast,  we  generated  MCF-7 
clones  with  reduced  HuR  levels  (93%)  using  lentiviral 
shRNA  (Figure  7B).  Western  blot  analysis  of  MCF-7 
cells  with  over-expression  of  HuR  reveals  modest 
increases  in  CD9  protein  levels  (Figure  7C).  There  are 
also  modest  decreases  in  CD9  protein  expression  in 
MCF-7  with  reduced  HuR  levels  (Figure  7C).  mRNA 
levels  of  CD9  and  CALM2  are  essentially  unchanged  in 
MCF-7  cells  which  over-express  HuR  (Figure  7D).  As 
expected,  HuR  knock-down  in  MCF-7  cells  using  lenti¬ 
viral  shRNA  resulted  in  significant  reductions  in  both 
CD9  and  CALM2  mRNA  levels  (Figure  7E).  The  right 
panel  in  Figure  7E  indicates  efficiency  of  HuR  mRNA 
knock-down  which  is  consistent  with  the  protein  data 
(Figure  7B).  These  results  are  summarized  in  Figure  7F. 
There  are  no  significant  changes  seen  in  CD9  mRNA 
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Figure  6  HuR  differentially  regulates  CD9  and  CALM2  in  MB-231.  A.  Epitope  HA  tagged  HuR  is  over-expressed  by  142%  and  138% 
respectively,  in  stably  transfected  clones  4E1  and  5F1,  as  compared  to  empty  vector  (EV)  control  clone  2C7.  B.  HuR  knock-down  using  lentiviral 
short  hairpin  (sh)  RNA  H760  results  in  a  94%  reduction  in  steady  state  levels  of  protein  in  clone  A7  (LL  =  lentilox  control).  C.  HuR  over-expression 
results  in  a  40%  reduction  in  CD9  protein  levels  as  assayed  by  Western  analysis;  however,  HuR  knock-down  using  lentiviral  shRNA  results  in  an 
increase  from  100%  to  228%  of  CD9  levels.  D.  Over-expression  of  HuR  decreases  CD9  mRNA  levels  but  not  CALM2  expression.  Analysis  of  steady 
state  CD9  and  CALM2  mRNA  levels  by  quantitative  RT-PCR  reveals  significant  decreases  in  CD9  mRNA  levels,  whereas  CALM2  levels  are 
unaffected.  Although  CALM2  expression  appears  greater,  the  change  is  not  significant.  E.  Knocking  down  HuR  levels  by  shRNA  in  MB-231  cells 
shows  significant  increases  in  CD9  and  CALM2  mRNA  levels  by  quantitative  RT-PCR.  Decreased  levels  of  HuR  mRNA  validate  HuR  shRNA  knock¬ 
down.  F.  Graph  showing  the  effects  of  HuR  on  the  expression  of  CD9  mRNA.  HuR  over-expression  results  in  decreases  in  both  mRNA  and 
protein  levels,  though  the  decreases  are  greater  in  RNA.  Whereas,  HuR  knock-down  by  shRNA  results  in  significant  increases  at  both  the  mRNA 
and  protein  levels,  with  greater  change  at  transcript  levels.  The  dashed  line  represents  levels  in  control  cells.  Error  bars  represent  SEM.  p  value  is 
<  0.005;  N.S.  =  not  statistically  significant;  *  =  statistically  significant.  All  experiments  were  done  in  triplicate  (n  =  3). 


and  CD9  protein  for  HuR  over-expression.  There  is  a 
more  pronounced  knock-down,  however,  in  CD9  mRNA 
in  MCF-7  cells  with  reduced  HuR  levels. 

The  results  of  HuR  shRNA  knock-down  in  MCF-7  cells 
were  as  expected,  but  opposite  of  those  seen  for  MB-231 
cells.  Steady-state  mRNA  levels  of  CD9  and  CALM2 
mRNAs  decreased,  consistent  with  the  hypothesis  that 
HuR  generally  stabilizes  its  mRNA  targets.  One  possible 
explanation  of  these  disparate  results  is  different  levels  of 


total  cellular  or  cytoplasmic  HuR.  We  performed  nuclear 
and  cytoplasmic  fractionation  (Additional  File  3,  Figure 
S3).  These  results  demonstrated  modest  (approximately 
10%)  greater  cytoplasmic  levels  of  HuR  in  MB-231  cells 
as  compared  to  MCF-7.  The  total  cellular  HuR  levels  are 
very  similar  for  both  MB-231  and  MCF-7  cells.  Taken 
together,  these  results  indicated  that  HuR  appeared  to 
differentially  regulate  the  same  mRNAs  in  a  manner 
dependent  upon  the  cellular  milieu. 
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Figure  7  Effects  of  over-expressing  or  reducing  HuR  on  CD9  and  CALM2  expression  in  MCF-7  cells.  A  Western  analysis  of  HuR  over¬ 
expression  in  heterogenous  population  of  cells  reveals  approximately  10%  over-expression.  B.  Lentiviral  HuR  shRNA  efficiently  knocks  down  HuR 
protein  by  over  90%.  C.  HuR  over-  expression  and  under-expression  results  in  small  changes  in  CD9  protein  levels  in  MCF-7  cells.  D.  Levels  of 
both  CD9  and  CALM2  mRNAs  are  unchanged  in  cells  which  over-express  HuR;  whereas  lentiviral  knock-down  of  HuR  in  MCF-7  cells  results  in 
decreases  in  steady-state  mRNA  levels  (E).  The  graph  in  (F)  shows  minimal  changes  in  CD9  mRNA  and  protein  levels  in  HuR  over-expressing 
MCF-7  cells.  The  CD9  mRNA  levels,  however,  are  more  affected  in  HuR  knock-down.  P  value  is  <  0.005;  N.S.  =  not  statistically  significant;  *  = 
statistically  significant. 


Discussion 

We  utilized  RIP-Chip  technologies  to  define  differentially 
regulated  HuR  genes  in  ER+  and  ER-  breast  cancer.  To 
our  knowledge,  this  is  the  first  report  of  a  side-by-side 
genome-wide  comparison  of  HuR-associated  targets  in 
wild  type  ER+  and  ER-  breast  cancer  cells.  Our  findings 
indicated  that  HuR  interacts  with  small  subsets  of  genes 
in  breast  cancer,  out  of  the  possible  8%  of  human  genes 
possessing  AREs  which  are  potential  HuR  targets.  Three 
broad  categories  of  HuR  targets  were  identified.  First, 


there  was  a  subset  of  targets  only  found  in  ER+  breast 
cancer.  Second,  there  was  a  unique  subset  of  HuR  targets 
found  only  in  ER-  breast  cancer.  A  third  subset  consisted 
of  HuR-associated  mRNAs  common  to  both  forms  of 
breast  cancer,  many  of  which  were  previously  described 
as  having  roles  in  cancer. 

We  selected  and  validated  two  HuR  targets,  CD9  and 
CALM2  mRNAs,  which  were  found  in  high  abundance  in 
both  types  of  breast  cancer.  Initially,  we  employed  the  pre¬ 
viously  developed  “heat  map”  signature  of  HuR  binding  to 
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gain  insight  into  putative  HuR  target  sequences  [30] .  HuR 
binding  was  verified  by  HuR  immunoprecipitations,  and 
analyzed  by  quantitative  RT-PCR  and  biotin  pull-downs. 
Both  CD9  and  CALM2  mRNAs  were  enriched  in  HuR 
RIPs  compared  to  isotype  control  IP  reactions.  Biotin  pull- 
downs  verified  the  binding  of  HuR  protein  specifically  to 
the  3’UTR  regions  of  both  mRNAs,  as  had  been  predicted. 
CD9  is  a  tetraspanin  molecule  which  plays  important  roles 
in  cellular  development,  activation,  growth  and  motility.  It 
has  been  implicated  in  a  variety  of  cancers,  including  but 
not  limited  to  gastric  cancers  and  B  cell  acute  leukemia 
[58-60]. 

The  role  of  CALM2  in  cancer  is  less  well  understood 
but  may  be  linked  to  cancer  since  it  is  involved  in  con¬ 
trolling  calcium  signaling  [61,62].  There  are  three  CAL¬ 
MODULIN  genes  ( CALM1 ,  CALM2  and  CALM3)  highly 
expressed  in  both  MB-231  and  MCF-7  cell  lines  (Addi¬ 
tional  File  4,  Figure  SI).  Interestingly,  although  they  are 
encoded  by  different  genes  at  different  chromosomal 
locations,  all  three  encode  the  same  open  reading  frame 
but  differ  in  the  5’  and  3’  UTRs  [[61,63],  and  [64]].  Only 
CALM2  mRNA  interacts  with  HuR  by  RIP  analysis. 
Moreover,  previously  published  reports  have  indicated 
the  necessity  of  knocking  down  all  three  CALMODULIN 
mRNAs  by  siRNA  to  achieve  knock-down  of  the  protein 
[61].  We  conclude  that  there  may  be  differential  HuR 
associated  regulation  of  these  CALMODULIN  genes  in 
breast  cancer,  even  though  the  mechanism  needs  to  be 
further  delineated. 

Surprisingly,  the  regulation  of  both  CD9  and  CALM2 
target  genes  was  dependent  upon  the  cellular  milieu.  To 
test  the  functional  consequences  of  HuR  binding  to 
these  two  transcripts,  we  prepared  cells  that  stably 
expressed  higher  or  lower  HuR,  compared  to  the  parent 
cells,  in  both  ER+  and  ER-  breast  cancer  cell  lines.  HuR 
appears  to  differentially  regulate  the  expression  of  CD9 
in  opposite  directions  in  the  two  different  forms  of 
breast  cancer.  Specifically,  HuR  over-expression  in  ER- 
breast  cancer  (MB-231)  paradoxically  decreased  CD9 
mRNA  and  protein  levels,  whereas  HuR  knock-down 
increased  the  CD9  mRNA  levels.  This  is  the  opposite  of 
what  is  predicted  for  most  HuR  targets,  since  HuR  is 
thought  to  stabilize  its  mRNA  targets  and  often 
increases  their  translation.  There  did  not  seem  to  be 
similar  effects  upon  CALM2  expression.  As  expected, 
knock-down  of  HuR  by  shRNA  decreased  expression  of 
CD9  and  CALM2  in  ER+  breast  cancer  (MCF-7). 
Though  there  are  differences  in  cytoplasmic  HuR  levels 
in  MB-231  cells  as  compared  with  MCF-7,  these  are 
modest  (10%).  This  is  in  keeping,  however,  with  obser¬ 
vations  that  MB-231  cells  are  more  undifferentiated  and 
more  aggressive. 

Moreover,  analysis  of  HuR-associated  mRNAs  in  both 
ER+  and  ER-  breast  cancer  revealed  three  broad 


categories  of  genes.  First,  there  were  well  known  cancer 
genes,  such  as  PTMA,  which  are  regulated  by  HuR  [27]. 
Second,  there  were  cancer-related  genes,  such  as  CD9 
and  CALMODULIN,  which  were  not  known  to  be  HuR 
regulated  until  this  report.  Third,  there  were  other 
genes  identified  by  HuR  association  with  unknown  can¬ 
cer  function.  These  could  potentially  represent  novel 
cancer  targets.  Additional  proof  of  HuR  involvement 
with  other  known  cancer  genes,  such  as  CD44  and 
GATA-3,  may  represent  novel  insights  into  the  mechan¬ 
isms  of  regulation  of  these  cancer  targets  (see  Additional 
Files).  These  results  may  therefore  advance  the  field  by 
shedding  insights  into  posttranscriptional  regulation  of 
known  and  perhaps  unknown  cancer  target  genes. 

Though  the  exact  mechanisms  of  HuR  differential  reg¬ 
ulation  of  CD9  and  CALM2  are  presently  unclear,  it 
may  involve  microRNA  (miRNA)  regulation.  In  a  recent 
report,  we  described  the  recruitment  by  HuR  of  miRNA 
let-7  to  translationally  silence  C-MYC  expression  [46].  It 
is  clear  from  the  findings  of  laboratories  headed  by  Fili- 
powicz,  Steitz  and  other  investigators,  that  RBPs  and 
miRNAs  are  involved  in  intricate  associations  to  affect 
downstream  translational  suppression  or  activation  of 
target  mRNAs  to  help  meet  cellular  needs  [65,66].  Sharp 
and  colleagues  proposed  that  different  interactions 
between  RBPs  and  miRNAs  may  have  evolved  as  a  pro¬ 
tective  mechanism  for  the  cell  against  environmental 
stress  [67]. 

A  remaining  question  is  why  HuR  selectively  binds  to 
certain  genes  containing  AREs.  Our  previous  work  has 
demonstrated  the  role  that  HuR  plays  in  myogenesis  by 
stabilizing  the  expression  of  three  critical  genes  involved 
in  myogenesis:  MYOD,  MYOGENIN,  and  p21cipl  [68], 
HuR  over-expression  results  in  precocious  muscle  differ¬ 
entiation  and  HuR  siRNA  knock-down  prevents  muscle 
differentiation  [69].  It  is  highly  probable  that  there  are 
more  than  three  HuR  targets  inside  these  cells.  A  speci¬ 
fic  phenotype  potentially  arises  when  HuR  levels  are 
altered  which  may  involve  interactions  with  miRNAs, 
although  this  theory  needs  to  be  fully  investigated. 

Our  findings  share  some  similarity  to  earlier  reports  of 
HuR  RIP-Chip  analysis  of  MCF-7  cells  stably  transfected 
with  MCT-1  [45].  These  analyses,  however,  were  not  gen¬ 
ome-wide  and  employed  transfected  cells.  Nevertheless, 
thrombospondin,  a  known  important  anti-angiogenic  fac¬ 
tor,  was  identified  as  a  HuR-regulated  target.  Combined 
with  earlier  reports  of  the  role  of  HuR  in  regulating, 
VEGF-a  and  HIFla,  HuR  may  be  controlling  a  “posttran¬ 
scriptional  mini-operon”  involved  in  angiogenesis 
[29,32,70].  Further  studies  are  being  conducted  in  our 
laboratory  to  investigate  the  role  of  HuR  in  breast  cancer 
angiogenesis  using  xenograft  animal  models.  It  will  be  par¬ 
ticularly  important  to  test  the  role  of  HuR  upon  CD9  and 
CALM2  expression  in  breast  tumors  in  vivo. 
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Posttranscriptional  gene  regulation  is  increasingly 
being  appreciated  as  a  driver  of  malignant  transforma¬ 
tion.  The  roles  of  both  RBPs  and  miRNAs  (so-called 
oncomirs)  are  being  recognized  in  cancer  [71].  Many 
reports  have  described  alterations  in  miRNA  expression 
profile  and  function  as  contributing  to  breast  cancer 
malignant  transformation  and  metastasis  [72-75].  HuR 
RIP-Chip  analysis  may  shed  further  light  into  malignant 
breast  cancer  transformation  by  identifying  HuR  asso¬ 
ciated  mRNAs. 

We  believe  that  there  are  potential  applications  for 
tamoxifen  resistance  as  well.  Keen  and  colleagues  have 
described  a  potential  mechanistic  link  between  HuR 
expression  and  tamoxifen  drug  resistance  [76].  As  breast 
cancer  cells  acquire  tamoxifen  resistance,  there  are 
increased  levels  of  cytoplasmic  HuR  expression. 
Increased  cytoplasmic  HuR  levels  have  previously  been 
described  in  situations  where  HuR  actively  influences 
expression  of  cytoplasmic  targets  [18,47,48].  Drug  resis¬ 
tance  could  be  reversed  by  using  siRNA  to  knock-down 
HuR  expression,  whereas  exogenous  over-expression  of 
HuR  could  cause  cells  to  become  resistant  to  tamoxifen. 
We  therefore  propose  that  HuR  may  be  coordinately 
regulating  genes  which  may  allow  a  cell  to  acquire 
tamoxifen  resistance.  It  will  be  interesting  to  further 
investigate  HuR-associated  target  genes  in  ER+  cells  in 
this  light. 

Conclusion 

In  summary,  using  RIP-Chip  analysis,  we  have  per¬ 
formed  for  the  first  time  a  genome-wide  comparison  of 
HuR-associated  targets  in  wild  type  ER+  and  ER-  breast 
cancer.  We  have  identified  novel  HuR  targets  and  have 
gained  insight  into  the  role  HuR  plays  in  regulating 
known  cancer  genes.  We  found  distinct,  differentially 
expressed  subsets  of  HuR  cancer  related  genes  in  ER+ 
and  ER-  breast  cancer  cell  lines.  Based  on  our  observa¬ 
tions,  the  enhanced  expression  of  these  mRNA  subsets 
by  HuR  can  influence  many  of  the  acquired  capabilities 
of  cancer  cells.  Further  investigation  into  HuR’s  role  in 
regulating  these  genes  may  provide  novel  insights  into 
breast  cancer  diagnosis  and  therapy. 

Appendix 

Microarray  Data  Preprocessing 

Data  quality  was  examined  by  looking  at  quality  controls 
metrics  produced  by  Illumina’s  software  (BeadStudio 
v3. 1.3.0,  Gene  Expression  Module  3.2.7).  The  data  were 
then  exported  for  further  analyses  in  R.  Image  plots  of 
each  array  were  examined  for  spatial  artifacts,  and  there 
was  no  evidence  of  systematic  effects  indicative  of  tech¬ 
nical  problems  with  the  arrays.  Within  limma,  quantile 
normalization  was  used  for  between  chip  normalization. 
Finally,  quality  control  statistics  were  computed  using  a 


variety  of  Illumina’s  internal  control  probes  that  are 
replicated  on  each  array.  Any  probes  which  were  consid¬ 
ered  “not  detectable”  across  all  samples  were  excluded 
from  further  statistical  analyses  in  order  to  reduce  false 
positives.  The  determination  of  “not  detectable”  was 
based  upon  the  BeadStudio  computed  detection  /7-value 
being  greater  than  1%. 

Gene  Ontology  Gene  Universe 

In  defining  the  gene  universe  for  the  analysis,  non-speci¬ 
fic  filtering  was  used  to  increase  statistical  power  with¬ 
out  biasing  the  results.  We  started  with  all  probes  on 
the  Illumina  array  which  had  both  an  Entrez  gene  iden¬ 
tifier  [77]  and  a  GO  annotation,  as  provided  in  the  lumi- 
HumanAll.db  [78]  annotation  data  package  and  GO.db 
[79]  annotation  maps  (built  using  data  obtained  from 
NCBI  on  4/2/08).  This  set  was  then  reduced  by  exclud¬ 
ing  probes  that  exhibited  little  variability  (interquartile 
range  (IQR)  of  <0.1  on  log2  scale)  across  all  samples 
because  such  probes  are  generally  not  informative. 
Finally,  for  probes  that  mapped  to  the  same  Entrez  iden¬ 
tifier,  a  single  probe  was  chosen  in  order  to  insure  a 
surjective  map  from  probe  IDs  to  GO  categories  (via 
Entrez  identifiers).  This  was  necessary  to  avoid  redun¬ 
dantly  counting  GO  categories  which  produces  false 
positives.  Probes  with  the  largest  IQR  were  chosen  to  be 
associated  with  an  Entrez  identifier. 


Additional  file  1:  Figure  S2.  Table  of  complete  GO  analysis.  Listing  of 
HuR-associated  genes  with  odds  ratios  and  functional  categories. 

Additional  file  2:  Table  SI.  HuR  targets  five  fold  or  greater.  Listing  of 
HuR-associated  mRNAs  in  MB-231  and  MCF-7  cell  lines. 

Additional  file  3:  Figure  S3.  Total  cellular  levels  of  HuR  are  similar 
in  MB-231  and  MCF-7  cells.  Nuclear  and  cytoplasmic  separation  was 
performed  to  measure  levels  of  HuR  in  different  compartments  of  MB- 
231  and  MCF-7  cells.  Total  cellular  HuR  levels  were  very  similar,  whereas 
there  was  a  small  (10%)  increase  in  HuR  cytoplasmic  levels  in  MB-231 
cells  as  compared  to  MCF-7.  Absence  of  p-tubulin  staining  demonstrates 
integrity  of  isolation  as  there  should  not  be  p-tubulin  in  the  nuclear 
fraction.  Bands  were  measured  by  densitometry  and  normalized  to  p- 
tubulin  controls.  (T  =  total  cellular  lysate;  C  =  cytoplasmic  lysate,  N  = 
nuclear  lysate). 

Additional  file  4:  Figure  SI.  Relative  baseline  values  of  CALM 7, 
CALM2,  CALM3,  and  CD-9  mRNAs  in  ER+  and  ER-  cells.  Quantitative 
RT-PCR  performed  on  mRNA  extracted  from  cell  lysates  showing  relative 
levels  of  CALM1,  CALM2,  CALM3 ,  and  CD-9  mRNAs  in  MB-231  and  MCF-7 
breast  cancer  cells.  All  values  were  normalized  to  GAPDH  mRNA.  All 
experiments  were  done  in  triplicate  (n  =  3)  except  for  CALM3  (n  =  2). 


Abbreviations 

(ER-):  Estrogen  receptor  negative;  (ER+):  estrogen  receptor  positive;  (RIP): 
RNA  immunoprecipitation;  (RIP-Chip):  RNA  immunoprecipitation  applied  to 
microarrays;  (3'  UTR):  3'  untranslated  region;  ELAV1 :  (embryonic  lethal 
abnormal  vision  1). 
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Interactions  between  RNA  binding  proteins  (RBPs)  and  genes  are  not  well  understood,  especially  in  regulation  of 
angiogenesis.  The  RBP  HuR  binds  to  the  AU-rich  (ARE)  regions  of  labile  mRNAs,  facilitating  their  translation  into  protein 
and  has  been  hypothesized  to  be  a  tumor-maintenance  gene.  Elevated  levels  of  cytoplasmic  HuR  directly  correlate  with 
increased  invasiveness  and  poor  prognosis  for  many  cancers,  including  those  of  the  breast.  HuR  controls  the  expression 
of  multiple  genes  involved  in  angiogenesis  including  VEGFa,  HIFIa  and  thrombospondin  1  (TSP1).  We  investigated  the 
role  of  HuR  in  estrogen  receptor  negative  (ER  )  breast  cancer.  MDA-MB-231  cells  with  higher  levels  of  HuR  have  alterations 
in  cell  cycle  kinetics  and  faster  growth.  Unexpectedly,  HuR  overexpression  significantly  interfered  with  tumor  growth  in 
orthotopic  mouse  models.  The  putative  mechanism  seems  to  be  an  anti-angiogenetic  effect  by  increasing  expression 
of  TSP1  but  also  surprisingly,  downregulating  VEGF,  a  target  which  HuR  normally  increases.  Our  findings  reveal  that  HuR 
may  be  regulating  a  cluster  of  genes  involved  in  blood  vessel  formation  which  controls  tumor  angiogenesis.  An  approach 
of  modulating  HuR  levels  may  overcome  limitations  associated  with  monotherapies  targeting  tumor  vessel  formation. 


Introduction 

Although  most  breast  cancers  are  estrogen  receptor  positive 
(ER^),  approximately  15%  are  estrogen  receptor  negative  (ER).1 
In  general,  ER'  tumors  are  more  aggressive,  have  a  poor  clinical 
outcome  and  dismal  survival  rates  that  disproportionately  affect 
lower  income  and  minority  women.  The  poor  prognosis  is  partly 
attributed  to  the  fact  that  these  women  are  limited  to  conven¬ 
tional  cytotoxic  chemotherapy.  This  is  in  distinct  contrast  to  ER+ 
tumors,  which  can  be  treated  with  hormonal  therapies  such  as 
tamoxifen  or  aromatase  inhibitors.  Moreover,  the  presence  of 
the  growth  factor  HER2  allows  oncologists  to  employ  targeted 
therapy  using  HER2  receptor-specific  drugs. 

The  role  of  posttranscriptional  gene  regulation  in  cancer  is 
now  widely  appreciated,  as  the  contributions  of  RNA  binding 
proteins  (RBPs)  and  microRNAs  (miRNAs)  are  becoming  more 
apparent.2'8  The  posttranscriptional  operon  hypothesis  states 
that  RBPs  and  miRNAs  are  coordinately  regulating  expres¬ 
sion  of  biologically  related  mRNAs.  This  has  generated  inter¬ 
est  in  identification  of  additional  gene  products  in  the  acquired 


capabilities  model  of  malignant  transformation.9  The  RBP  HuR, 
a  paraneoplastic  antigen,  overexpressed  in  many  malignancies 
including  breast  cancer,  has  been  implicated  as  an  important 
RBP  which  may  function  as  a  tumor  maintenance  gene,  facilitat¬ 
ing  malignant  transformation.10'17  HuR  has  been  shown  to  regu¬ 
late  genes  in  multiple  areas  of  the  acquired  capabilities  model, 
including  two  pivotal  genes  involved  in  angiogenesis,  VEGF  and 
HIFlct. 18-22  Increased  HuR  cytoplasmic  expression  directly  cor¬ 
relates  with  severity  and  aggressiveness  of  many  cancers,  includ¬ 
ing  those  of  the  breast.23'25  HuR  has  been  demonstrated  to  be  an 
important  prognostic  factor  in  familial  breast  cancer  patients.23  24 
Recently,  investigators  altered  the  subcellular  distribution  of 
HuR  and  documented  a  decrease  in  mRNA  stability,  subsequent 
ER  levels  and  susceptibility  to  tamoxifen.26 

In  this  study  we  investigated  the  role  of  HuR  in  ER'  breast 
cancer  by  gain-of-function  ER'  (MDA-MB-231)  cell  lines  which 
overexpressed  epitope  tagged  hemagglutinin  (HA)  HuR.  As 
expected  and  consistent  with  published  results  for  other  can¬ 
cer  types,  HuR  overexpression  in  vitro  resulted  in  increasing 
MDA-MB-231  cellular  growth  rates  and  alterations  in  cell  cycle 
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Figure  1.  Overexpression  of  HA-HuR  in  MDA-MB-231  cancer  cells  increases  cell  growth  and  alters  cell  cycle  kinetics  in  vitro  while  inhibiting  tumor 
growth  in  vivo.  (A)  MDA-MB-231  cells  transfected  with  the  pZeoSV2  vector  expressing  HA-HuR,  selected  with  Zeocin  and  cloned  by  limiting  dilution 
express  HA-HuR  compared  to  a  pZeoSV2  empty  vector  control.  MDA-MB-231  clones  4E1  and  5F1  expressed  42%  and  38%  more  HuR  than  control  as 
shown  in  representative  western  blot.  Bar  graph  represents  mean  HuR  percent  overexpression  from  three  different  western  blots.  (B)  Both  clones 
expressing  HA-HuR  proliferated  significantly  faster  than  the  empty  vector  control  in  vitro.  (C)  Overexpression  of  HA-HuR  increased  cells  in  G^G,  cell 
cycle,  from  57.70%  to  67.67%.  Overexpression  of  HA-HuR  decreased  cells  in  G2/M  phase  by  27.19%  to  18.29%  but  did  not  significantly  alter  cells  in  S 
phase  as  measured  by  DNA  content.  (D)  MDA-MB-231  HA-HuR  4E1  showed  significantly  reduced  tumor  volume  (mm3)  and  growth  starting  at  two 
weeks  post-inoculation  and  continuing  for  five  weeks  when  compared  to  empty  vector  MDA-MB-231  controls  as  measured  by  both  MRI  and  calipers. 
For  tumor  experiments  nine  animals  per  group  were  used.  For  both  counting  assay  and  tumor  volumes  data  represent  mean  value  ±  SEM.  p  <  0.05. 
Flow  cytometry  data  represent  mean  value  ±  SEM  from  n  =  4  separate  experiments  done  in  triplicate,  p  <  0.05. 


kinetics.  Surprisingly,  when  these  cells  were  used  in  vivo  as  xeno¬ 
grafts  in  athymic  mice,  a  90%  reduction  in  growth  rate  was  seen, 
as  compared  to  control  groups.  The  presumptive  mechanism 
was  an  anti-angiogenic  affect  by  increasing  the  expression  of 
thrombospondin  1  (TSP1),  a  gene  known  to  be  anti-angiogenic. 
However,  unexpectedly,  there  were  also  decreases  in  VEGF 
expression.  Thus,  our  results  indicate  that  HuR  overexpression 
in  ER  breast  cancer  can  potently  inhibit  tumor  growth  by  poten¬ 
tially  blocking  angiogenesis. 

Results 

Overexpression  of  HuR  in  MDA-MB-231  cells  increases 
growth  rates  and  alters  cell  cycle  kinetics.  To  study  the  role  of 
HuR  expression  in  MDA-MB-231  ER  breast  cancer  we  made 
individual  clones  which  overexpressed  either  epitope  tagged  (HA) 
HuR  or  empty  vector  (EV)  control  and  measured  growth  rates 
and  cell  cycle  kinetics  (Fig.  1A).  Overexpression  of  HA-HuR 


resulted  in  accelerated  cellular  growth  as  determined  by  counting 
(Fig.  IB).  When  the  cells  were  stained  with  propidium  iodide, 
we  noted  an  alteration  in  cell  cycle  kinetics  (Fig.  1C).  HA-HuR 
overexpressing  cells  had  increased  amounts  of  cells  in  Gt  (67 
vs.  57%),  as  compared  with  empty  vector  controls.  Conversely, 
HuR  overexpression  also  resulted  in  a  compensatory  decrease  in 
G2/M  percentages  (18%  vs.  27%).  We  concluded,  as  expected, 
that  HuR  overexpression  resulted  in  increases  in  growth  rates  of 
MDA-MB-231  cells.  We  then  investigated  the  effects  of  HuR 
overexpression  in  vivo,  using  orthotopic  xenograft  animal  models. 

HuR  overexpression  results  in  significantly  reduced  tumor 
growth  and  mass.  The  clones  used  in  Figure  1,  empty  vector 
(2C7)  and  overexpresser  HA-HuR  (4E1),  were  injected  into  the 
contralateral  mammary  fat  pads  of  athymic  nude  mice.  Tumor 
growth  was  assessed  weekly  by  caliper  measurements  and  fol¬ 
lowed  in  vivo  by  MRI  scan.  Surprisingly,  tumors  overexpressing 
HuR  had  significant  inhibition  of  growth,  whereas  EV  control 
tumors  increased  significantly  in  tumor  volume  (Fig.  ID).  The 
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Figure  2.  Overexpression  of  HA-HuR  in  MDA-MB-231  cancer  cells  inhibits  tumor  growth  in  athymic  nude  mice.  (A)  Repeat  experiments  comparing 
MDA-MB-231  HA-HuR  4E1  with  both  wild-type  MDA-MB-231  and  empty  vector  MDA-MB-231  confirmed  HuR  overexpression  reduced  tumor  volume 
(mm3)  and  growth  starting  at  two  weeks  post-inoculation  and  continuing  for  five  weeks  as  measured  by  calipers.  (B)  Tumors  overexpressing  HA-HuR 
had  significantly  less  mass  after  harvest  42  days  post-inoculation  when  compared  to  the  WT  or  empty  vector  (EV)  controls.  (C)  MRI  comparing  largest 
sections  for  each  tumor  showed  significantly  smaller  tumors  in  the  HA-HuR  overexpressing  tumors  when  compared  to  EV  control  tumors.  (D)  Repre¬ 
sentative  cross  sections  of  tumors  showed  that  those  formed  by  inoculation  with  HA-HuR  resembled  a  gelatin-like  capsule,  were  significantly  smaller 
and  more  homogeneous  than  those  formed  by  inoculation  with  empty  vector.  Hematoxylin  and  eosin  stain  revealed  poorly  differentiated  carcinomas 
with  similar  morphology  and  lack  of  inflammatory  cells  in  both  HA-HuR  tumors  and  EV  control  tumors.  Five  animals  were  used  per  group  in  HA-HuR, 
empty  vector  and  wild-type  control  groups.  Experiments  were  repeated  with  similar  results  using  a  different  clone,  5F1,  (see  Supporting  Information 
Fig.  S3).  Data  represent  mean  value  ±  SEM.  p  <  0.05;  in  photomicrographs  bar  =  27  microns. 


tumors  were  removed  from  the  animals  on  day  42.  Histological 
staining,  western  blotting  and  RT-PCR  were  performed  to  deter¬ 
mine  whether  any  cells  remained.  As  seen  in  Supplementary 
Data,  both  control  (2C7)  and  HuR  (4E1)  tumors  had  intact 
human  GAPDH  mRNA  and  tumor  4E1  expressed  the  HA-HuR 
transgene.  Furthermore,  the  HA-HuR  tumors  harvested  on  day 
42  still  expressed  HA-HuR  protein  (Fig.  SI).  We  reestablished 
cell  lines  from  tumors  removed  from  animals  to  study  their 
growth  rates.  We  found  that  these  reestablished  cells  had  similar 
growth  rates  to  parental  cells  prior  to  transplantation  (Fig.  S2). 

The  experiments  were  repeated  with  wild-type,  parental 
MDA-MB-231  cells  to  validate  these  findings.  As  seen  in  Figure 
2,  both  parental  and  control  cells  (2C7)  grew  at  similar  rates, 
whereas  tumors  which  overexpressed  HuR  (4E1)  did  not  appear 
to  grow  (Fig.  2A).  Tumor  mass  was  assessed  and  showed  that 
HuR  overexpression  resulted  in  a  90%  reduction  in  growth 
(Fig.  2B).  These  results  were  confirmed  by  MRI  scans,  gross 
photographs  and  microscopy  (Fig.  2C  and  D).  HA-HuR  tumors 
appeared  to  be  a  gelatin-like  capsule  and  the  control  tumors 
were  a  solid  round  mass.  Cross  sections  of  both  revealed  that  the 
HA-HuR  tumors  had  a  smooth,  homogeneous  and  glistening 
surface,  whereas  the  control  tumors  had  a  heterogeneous,  yel¬ 
low-white  surface  with  a  necrotic  center  (Fig.  2D).  Both  tumors 


contained  viable  cancer  cells  and  similar  morphology  determined 
to  be  moderately  to  poorly  differentiated  carcinoma,  consistent 
with  the  implanted  MDA-MB-231  cells  (Fig.  2D).  We  con¬ 
cluded  that  HuR  overexpression  resulted  in  significantly  smaller 
ER'  tumors  in  animals. 

To  verify  that  these  results  were  not  clonal,  we  repeated  the 
orthotopic  tumor  injection  experiments  using  a  second  HuR 
overexpressing  clone  (5F1)  which  in  vitro  grew  similarly  to  the 
original  overexpressing  clone,  4E1  (Fig.  IB).  As  seen  in  Figure 

53,  tumor  5F1  also  exhibited  retarded  growth  rates  as  com¬ 
pared  to  empty  vector  controls  and  had  a  90%  reduction  in  mass 
(Fig.  S3).  Taken  together,  these  findings  demonstrated  that  HuR 
overexpression  in  MDA-MB-231  cells  resulted  in  significant 
reductions  in  tumor  growth  in  a  clonal  independent  fashion. 

Gene  ontology  (GO)  analysis  of  overexpressed  genes  in 
HuR  overexpressing  cells.  In  order  to  better  understand  the 
genes  which  may  be  involved  in  altering  tumor  growth  in  HuR 
overexpressing  MDA-MB-231  cells,  we  performed  genome  wide 
microarray  analysis  (Fig.  S4  and  Table  1).  As  seen  in  Figure 

54,  many  genes  were  over-represented  (reflected  in  large  odds 
ratios)  that  deal  with  both  biological  processes  as  well  as  molecu¬ 
lar  function.  Given  the  large  numbers  of  genes,  we  decided  to 
investigate  the  following  three  potential  mechanisms  to  explain 
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Table  1.  Genes  upregulated  in  tumors  overexpressing  HuR 


Fold  change 

Entrez  ID 

Gene  name 

21.52 

79191 

IRX3/iroquois  homeobox  3 

16.87 

4256 

MGP/matrix  Gla  protein 

12.21 

1278 

COL1A2/collagen,  type  1,  alpha2 

9.35 

8714 

ABCC3/ATP-binding  cassette,  sub-family  C  (CFTR/MRP),  member  3 

7.58 

221476 

PI16/peptidase  inhibitor  16 

7.46 

972 

CD74/CD74  molecule,  major  histocompatibility  complex,  class  II  invariant  chain 

6.90 

2202 

EFEMPI/EGF-containing  fibulin-like  extracellular  matrix  protein  1 

6.84 

5265 

SERPINAI/serpin  peptidase  inhibitor,  clade  A  (alpha-1  antiproteinase,  antitrypsin),  member  1 

6.77 

3399 

ID3/inhibitor  of  DNA  binding  3,  dominant  negative  helix-loop-helix  protein 

6.63 

6275 

S100A4/S100  calcium  binding  protein  A4 

5.85 

57537 

SORCS2/sortilin-related  VPS10  domain  containing  receptor  2 

5.07 

6275 

S100A4 

5.07 

1831 

TSC22D3/TSC22  domain  family,  member  3 

4.87 

51313 

C4orf18/chromosome  4  open  reading  frame  18 

4.87 

4061 

LY6E/lymphocyte  antigen  6  complex,  locus  E 

4.77 

3108 

FILA-DMA/major  histocompatibility  complex,  class  II,  DM  alpha 

4.61 

10507 

SEMA4D/sema  domain,  immunoglobulin  domain  (Ig),  transmembrane  domain  (TM)  and 

short  cytoplasmic  domain,  (semaphorin)  4D 

4.46 

7052 

TGM2/transglutaminase  2  (C  polypeptide,  protein-glutamine-gamma-glutamyltransferase) 

4.46 

1728 

NQ01/NAD(P)Fi  dehydrogenase,  quinone  1 

4.43 

9022 

CLIC3/chloride  intracellular  channel  3 

4.42 

7057 

TFIBSI/thrombospondin  1 

4.36 

7857 

SCG2/secretogranin  II  (chromogranin  C) 

4.29 

4599 

MXI/myxovirus  (influenza  virus)  resistance  1,  interferon-inducible  protein  p78  (mouse) 

4.24 

22998 

LIMCH1/LIM  and  calponin  homology  domains  1 

4.14 

3371 

TNC/tenascin  C 

4.13 

6376 

CX3CL1/chemokine  (C-X3-C  motif)  ligand  1 

4.1 

199 

AIFI/allograft  inflammatory  factor  1 

4.05 

652 

BMP4/bone  morphogenetic  protein  4 

3.94 

7058 

TFIBS2/thrombospondin  2 

3.86 

26227 

PFIGDFI/phosphoglycerate  dehydrogenase 

3.78 

144165 

PRICKLEI/prickle  homolog  1  (Drosophila) 

3.76 

4147 

MATN2/matrilin  2 

3.75 

8140 

SLC7A5/solute  carrier  family  7  (cationic  amino  acid  transporter,  y+  system),  member  5 

3.54 

8460 

TPSTI/tyrosylprotein  sulfotransferase  1 

3.48 

64764 

CREB3L2/cAMP  responsive  element  binding  protein  3-like  2 

3.44 

3397 

IDI/inhibitor  of  DNA  binding  1,  dominant  negative  helix-loop-helix  protein 

3.42 

9516 

LITAF/lipopolysaccharide-induced  TNF  factor 

3.4 

55971 

BAIAP2L1/BAI1 -associated  protein  2-like  1 

3.36 

80063 

ATF7IP2/activating  transcription  factor  7  interacting  protein  2 

3.26 

1021 

CDK6/cyclin-dependent  kinase  6 

3.25 

7100 

TLR5/toll-like  receptor  5 

3.17 

94121 

SYTL4/synaptotagmin-like  4 

3.16 

8622 

PDE8B/phosphodiesterase  8B 

3.15 

55616 

DDEFLI/ArfGAP  with  SH3  domain,  ankyrin  repeat  and  PFH  domain  3 

Tumor  microarray  reveals  47  annotated  genes  up  regulated  in  the  HA-HuR  overexpressing  tumors  as  compared  to  EV  control  tumors.  47  genes  are  up 
regulated  3-fold  or  greater  in  the  HA-HuR  tumors  as  compared  to  EV  control  tumors  (false  discovery  rate  <1%),  and  also  have  a  probability  of  differen¬ 
tial  expression  >80%  based  on  a  Bayesian  analysis. 
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Table  1.  Genes  upregulated  in  tumors  overexpressing  HuR  (continued) 


3.13 

30818 

KCNIP3/Kv  channel  interacting  protein  3,  calsenilin 

3.13 

3430 

IFI35/interferon-induced  protein  35 

3.07 

29887 

SNXIO/sorting  nexin  10 

3.01 

23052 

ENDODI/endonuclease  domain  containing  1 

Tumor  microarray  reveals  47  annotated  genes  up  regulated  in  the  HA-HuR  overexpressing  tumors  as  compared  to  EV  control  tumors.  47  genes  are  up 
regulated  3-fold  or  greater  in  the  HA-HuR  tumors  as  compared  to  EV  control  tumors  (false  discovery  rate  <1%),  and  also  have  a  probability  of  differen¬ 
tial  expression  >80%  based  on  a  Bayesian  analysis. 


the  large  discrepancy  seen  in  tumor  growth:  (1)  increased 
apoptosis,  (2)  increases  in  senescence  and  (3)  alteration  in 
angiogenesis.  The  microarray  profiling,  however,  did  not  dis¬ 
tinguish  between  direct  and  indirect  HuR  targets.  Therefore, 
we  decided  to  further  investigate  known  HuR  targets  related  to 
cancer. 

Tumors  which  overexpress  HuR  have  decreased  angiogene¬ 
sis.  We  conducted  experiments  which  targeted  well  known  HuR 
target  genes  involved  in  angiogenesis  which  had  previously  been 
described  in  the  literature:  thrombospondin  1  (TSP1),  VEGF 
and  HIF  la. 21,22,27  We  therefore  performed  real-time  PCR  to  mea¬ 
sure  mRNA levels  ofTSPl,  VEGF  and  HIFla.  As  seen  in  Figure 
3A,  HuR  overexpression  caused  an  increase  in  TSP1  mRNA 
and  protein  (Fig.  3B)  but  surprisingly  decreases  in  steady-state 
VEGF  mRNA  and  protein  levels;  HIFla  steady-state  mRNA 
levels  appeared  unchanged  (Fig.  3A-C).  We  further  measured 
HIFla  protein  by  western  and  immunohistochemical  stain¬ 
ing,  which  were  unchanged  between  HuR  overexpressing  and 
control  tumors  (data  not  shown).  Furthermore,  RNA  immuno- 
precipitation  (RIP)  using  an  antibody  specific  to  HuR  and  an 
IgGl  antibody  control  revealed  that  HuR  interacts  with  TSP1 
and  VEGF  mRNA,  as  these  transcripts  were  38-  and  54-fold 
enriched,  respectively,  in  the  HuR  IP  pellet,  as  compared  with 
istotype  control.  TSP1  and  VEGF  mRNA  were  upregulated  62- 
and  59-fold  in  the  HuR  IP  from  the  HuR  overexpression  cells. 
GAPDH,  which  does  not  interact  with  HuR  showed  no  enrich¬ 
ment  in  either  HuR  IP  when  compared  to  the  IgGl  IP  (Fig.  4A). 
After  confirming  HuR  interacts  with  TSP1  and  VEGF  mRNA 
we  sought  to  gain  a  better  understanding  of  the  mechanism  by 
which  the  overexpression  of  HuR  leads  to  increases  in  TSP1 
mRNA  and  protein  and  decreases  in  VEGF  mRNA  and  protein. 
To  determine  whether  alterations  in  HuR  levels  were  altering 
TSP1  and  VEGF  mRNA  stability  we  performed  mRNA  decay 
assays  using  actinomycin  D.  The  VEGF  mRNA  stability  was  not 
altered  when  comparing  HuR  overexpressing  cells  with  empty 
vector  control  (Fig.  4B).  TSP1  mRNA  half-life  was  increased  in 
the  HuR  overexpressing  cells  when  compared  to  the  empty  vec¬ 
tor  control  cells  (Fig.  4C). 

We  next  examined  whether  alterations  in  apoptosis  could 
account  for  differences  in  tumor  growth.  To  determine  whether 
HuR  overexpression  was  altering  apoptosis  in  vitro  we  performed 
both  annexin  V  and  7-AAD  staining.  No  differences  in  amount 
of  cells  undergoing  apoptosis  between  the  HuR  overexpress- 
ers  and  EV  controls  were  seen  (Figs.  5A  and  S5).  In  vivo,  there 
were  no  increases  in  apoptosis  in  HuR  overexpression  tumors, 
as  compared  with  EV  control  tumors  harvested  on  day  42 


post-inoculation  (Figs.  5B  and  S6).  Increased  apoptosis  in  EV 
tumors  was  mostly  found  in  the  necrotic  centers.  To  account 
for  the  possibility  that  apoptosis  was  occurring  early  in  tumor 
formation,  we  repeated  the  experiments  and  harvested  tumors 
from  mice  14  days  post-inoculation.  No  increases  in  apopto¬ 
sis  were  seen  in  HA-HuR  tumors  as  compared  to  EV  control 
tumors  (Figs.  5B  and  S6).  The  tumors  from  mice  14  days  post¬ 
inoculation  were  only  beginning  to  be  visible.  Hematoxylin  and 
eosin  staining  revealed  tumor  morphology  consistent  with  breast 
adenocarcinoma  (Fig.  S7).  We  further  quantitated  the  amount 
of  blood  vessels  and  found  a  significant  decrease  in  vessel  forma¬ 
tion  in  HA-HuR  tumors  as  compared  with  EV  controls  (Fig.  5B 
and  C).  No  tortuous  vessels  were  seen.  These  non-functional 
vessels  are  sometimes  seen  when  there  are  perturbations  in 
the  DLL4-Notch  signaling  pathways.  Measurement  of  senes¬ 
cence  using  (3-galactosidase  staining  did  not  reveal  any  sig¬ 
nificant  differences  (Fig.  S8).  Taken  together,  we  concluded 
that  HuR  overexpression  in  MDA-MB-231  tumors  does  not 
appear  to  alter  cellular  apoptosis  but  may  interfere  with  angio¬ 
genesis  by  overexpression  of  an  anti-angiogenetic  factor,  TSP1 
and  by  downregulation  of  a  known  pro-angiogenic  factor, 
VEGF. 

Discussion 

In  ER'  cells  overexpressing  HuR,  there  are  increases  in  growth 
rates  in  vitro  as  well  as  alterations  in  cell  cycle  kinetics.  Specifically, 
MDA-MB-231  cells  overexpressing  HuR  have  increases  in  the  Gj 
phase  of  the  cell  cycle,  consistent  with  published  results.25  A  plau¬ 
sible  explanation  is  HuR  induced  stabilization  of  cyclin  Bl,  the 
pivotal  cyclin  involved  in  G2/M  transition.28  Surprisingly,  when 
these  same  cells  were  transplanted  into  athymic  nude  mice,  HuR 
overexpression  resulted  in  a  90%  growth  reduction  in  tumor 
size.  These  results  were  confirmed  by  measuring  tumor  volume 
and  mass.  Additionally,  the  results  were  validated  by  serial  MRI 
scans  during  the  experimental  period  and  further  confirmed  in 
two  independent  clones.  Empty  vector  (EV)  control  and  parental 
wild  type  MDA-MB-231  cells  grew  similarly  and  resulted  in  much 
larger  tumors  than  those  formed  by  HA-HuR  overexpressing  cells. 
Histological  staining  of  cross  sections  from  EV  and  HA-HuR 
tumors  revealed  them  to  be  poorly  differentiated  carcinomas,  con¬ 
sistent  with  ER'  breast  cancer.  Moreover,  the  HuR  transgene  was 
still  expressed  in  the  smaller  tumors  at  the  time  of  harvest  (day  42). 

We  searched  for  putative  mechanisms  to  explain  these  find¬ 
ings,  since  in  other  systems  HuR  overexpression  results  in  larger, 
more  robust  tumor  growth.25  No  inflammatory  infiltrates  were 


www.landesbioscience.com 


Cell  Cycle 


3341 


0) 

W  s 

c 

«  4 
£ 

o  , 

T3 


* 


I  Empty  Vector  control 


TSP1  VEGFa  HIFIa 

cell  tumor  tumor 

culture  *p<0.05 


g  HA  HuR  Empty  Vector 

Tumor  Control  Tumor 


TSP1 


p-Tubulin 

Percent  of 
EV  control 


176%  100% 


VEGFa 


P-Tubulin 

Percent  of 
EV  control 


HA  HuR 
Tumor 


67% 


Empty  Vector 
Control  Tumor 


100% 


Figure  3.  TSP1  is  upregulated  in  HA-HuR  tumor  and  VEGFa  is  down- 
regulated.  (A)  Real-time  PCR  indicates  TSP1  is  upregulated  in  tumors 
(5.44  fold)  and  cells  in  culture  (4.88  fold)  overexpressing  HA-HuR  when 
compared  to  EV  control  tumors  and  cells,  consistent  with  the  microar¬ 
ray  data.  VEGF  is  downregulated  (2.6  fold)  in  tumors  overexpressing 
ElA-EluR  when  compared  to  EV  controls.  HIFIa  mRNA  levels  did  not 
appreciably  change.  Change  in  gene  expression  was  determined  using 
the  comparative  CT  method  and  is  represented  as  fold  change  in  HA- 
HuR  tumors  as  compared  to  empty  vector  controls.  GAPDH  was  used 
as  an  endogenous  control.  (B)  Western  blot  for  TSP1  shows  increased 
protein  expression  of  TSP1  (76%)  in  the  HA-HuR  overexpressing  tumors 
when  compared  to  EV  control  tumors.  (C)  Western  blot  for  VEGF  shows 
decrease  protein  expression  by  23%  in  the  HA-HuR  overexpressing 
tumors  when  compared  to  EV  control  tumors  (representative  of  two 
independent  sets  of  tumors).  Data  represent  mean  value  ±  SEM  from 
n  =  3  separate  mice  done  in  triplicate,  p  <  0.05. 


seen  in  either  tumor.  As  expected,  there  was  increased  apopto¬ 
sis  in  the  centers  of  the  EV  control  tumors,  since  these  regions 
are  relatively  more  hypoxic.  Apoptosis  early  in  tumor  develop¬ 
ment  is  an  unlikely  mechanism  of  inhibited  tumor  growth  in 
the  tumors  overexpressing  HuR  as  we  did  not  see  any  altera¬ 
tions  in  apoptosis  14  days  post-inoculation  using  two  methods 
of  apoptosis  detection.  Furthermore,  significant  apoptosis  was 
not  observed  in  either  cell  lines  (in  vitro)  or  tumors  (in  vivo) 


overexpressing  HuR.  These  results  more  than  likely  rule  out 
apoptosis  as  a  mechanism  of  reduced  tumor  growth  in  HA-HuR 
tumors,  although  we  cannot  definitely  exclude  apoptosis  earlier 
than  14  days.  Since  the  tumors  continue  to  overexpress  HuR 
even  at  day  42,  it  would  seem  unlikely  that  HuR  overexpression 
would  be  the  etiology  of  cell  death.  We  also  ruled  out  senescence 
as  a  mechanism,  as  there  was  less  (3-galactosidase  staining  in  the 
smaller  tumors. 

Due  to  previously  published  data  on  the  role  of  HuR  in  con¬ 
trolling  angiogenesis  via  interactions  with  VEGF  and  HIFlot 
mRNAs,  we  investigated  the  relationship  between  HuR  over¬ 
expression  and  these  pro-angiogenetic  factors.21,22  There  was  a 
statistically  significant  decrease  in  VEGF  mRNA  and  protein 
expression  and  no  increases  in  HIFIa  mRNA  or  protein  expres¬ 
sion.  As  expected,  increased  HuR  expression  correlated  with 
increased  TSP1  mRNA  expression  and  protein  levels.  TSP1  is 
a  well  known  anti-angiogenetic  factor  and  has  been  described  to 
be  regulated  by  HuR.27  Enumeration  of  angiogenesis  by  CD34 
staining  confirmed  that  HuR  overexpression  resulted  in  signifi¬ 
cant  decreases  in  new  blood  vessel  formation.  This  decrease  in 
neovascularization  may,  in  part,  explain  why  tumors  formed  by 
HA-HuR  overexpression  were  much  smaller  than  EV  controls. 

Other  groups  have  shown  that  decreases  in  VEGF  expression 
in  a  variety  of  cancer  models  can  inhibit  angiogenesis  and  tumor 
growth.29'32  It  has  also  been  demonstrated  that  TSP1  can  mediate 
tumor  growth  suppression  via  blocking  angiogenesis,  consistent 
with  our  results.33,34 

As  previously  mentioned  HuR  has  been  described  to  stabilize 
TSP1  and  VEGF  mRNA  resulting  in  greater  levels  and  increased 
protein  expression.21,27  In  contradistinction  to  published  reports 
of  HuR  stabilization  of  VEGF  mRNA,  our  findings  indicate 
that  under  conditions  of  HuR  overexpression,  there  are  decreases 
in  both  VEGF  mRNA  and  protein  levels.  The  reasons  for  this 
discrepancy  are  presently  unknown,  though  the  results  may  be 
dependent  upon  levels  of  HuR  overexpression  and  hypoxia.  HuR’s 
relationship,  however,  with  HIFIa  is  more  complex.  HuR  binds 
to  AU-rich  regions  in  HIFIa  5'UTR,  rather  than  the  3'UTR 
(even  though  both  regions  of  the  molecule  possess  AREs),  and 
causes  a  translational  upregulation  in  HIFIa  protein  synthesis 
without  altering  mRNA  levels.  It  appears  HuR  overexpression  is 
not  affecting  HIFIa  protein  production  in  our  system.  HIFIa 
is  the  major  transcriptional  factor  for  VEGF  mRNA  transcrip¬ 
tion,  and  as  such,  this  is  an  area  of  intensive  investigation  in  our 
laboratory. 

The  exact  mechanisms  of  HuR-induced  anti-angiogenetic 
effects  are  not  completely  understood  but  may  involve  interac¬ 
tions  between  HuR  and  microRNAs.  In  a  recently  published 
report,  HuR  actively  recruited  let-7  miRNA  to  c-myc  mRNA, 
which  resulted  in  decreases  of  its  translation.35  Furthermore,  we 
do  not  have  any  direct  evidence  that  the  HA  tag  (located  at  the 
amino  terminal)  affects  the  distribution  or  targeting  of  HuR  to 
its  mRNA  targets.  Published  results  from  the  Kontoyonnis  lab 
used  the  same  epitope  tag  to  make  a  transgenic  mouse  overex¬ 
pressing  HA-HuR  in  macrophages  and  although  this  is  a  differ¬ 
ent  system  than  ours,  they  did  not  see  any  alterations  in  nuclear 
vs.  cytoplasmic  HuR  distribution  nor  binding  to  its  mRNA 
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Figure  4.  HuR  interacts  with  both  TSP1  and  VEGF  mRNAs  in  cells  overexpressing  HuR.  (A)  RNA  immuno- 
precipiation  indicates  both  TSP1  and  VEGF  mRNA  are  increased  in  the  HuR  IP  when  compared  to  IgGI 
control  IP  in  both  HA-HuR  overexpressing  cells  and  EV  control  cells.  (B)  Actinomycin  D  mRNA  stability 
assay  shows  VEGF  mRNA  half-life  was  not  altered  between  cells  overexpressing  HA-HuR  and  EV  control 
cells.  (C)  TSP1  mRNA  from  cells  overexpressing  HA-HuR  has  a  longer  half-life  than  TSP1  mRNA  from  EV 
control  cells.  For  RNA  immunoprecipitation,  data  represents  mean  value  ±  SEM.  p  <  0.05. 


targets.36  Moreover,  the  metastatic  potential  of  anti-angiogenetic 
MDA-MB-231  HuR  altered  cell  lines  needs  to  be  investigated, 
since  this  phenotype  results  in  patient  mortality  and  the  corre¬ 
lation  between  tumor  size  and  metastatic  potential  may  not  be 
direct.  We  have  identified  putative  markers  of  metastasis  which 
have  been  altered  in  vivo  in  these  cells  and  this  area  is  currently 
under  investigation. 

Interactions  between  RBPs  such  as  HuR  and  their  downstream 
mRNA  targets  are  complex.  Even  though  HuR  can  theoretically 
bind  8%  of  human  mRNAs  which  contain  AREs,  it  does  not 
necessarily  do  so.37'39  Published  reports  from  our  lab  and  others 
demonstrate  that  HuR  critically  regulates  myogenesis  by  control¬ 
ling  the  expression  of  three  central  genes  in  this  process.40,41  HuR 
siRNA  knock  down  prevents  muscle  formation,  whereas  HuR 
overexpression  results  in  precocious  differentiation.40,41  Yet,  it  is 
highly  likely  HuR  is  interacting  with  many  more  targets  inside 
these  cells.  The  exact  explanations  for  this  discrete  phenotype  are 
not  presently  understood.  Our  work  involving  targets  discovered 
by  RIP-Chip  in  cell  lines  under  normoxic  conditions  indicates 
that  HuR  may  indeed  regulate  the  same  targets  in  opposite  fash¬ 
ion,  dependent  upon  cellular  milieu,  though  the  reasons  are  not 
presently  clear.42  There  are  obvious  differences  when  performing 
these  experiments  under  the  hypoxic  conditions  which  operate  in 
animals  in  which  tumors  form.  The  published  data  has  both  indi¬ 
cated  that  HuR  overexpression  correlates  with  increased  as  well 
as  decreased  aggressiveness  in  breast  cancer.23,24,43  These  results 
would  lend  credence  to  the  idea  that  HuR  may  be  directing 


its  mRNA  targets  in  cell-specific 
fashion. 

In  summary,  our  findings  indi¬ 
cate  that  HuR  overexpression  in 
ER'  breast  cancer  results  in  a  con¬ 
comitant  increase  in  TSP1  and  a 
decrease  in  VEGF  expression  as  well 
as  substantial  decreases  in  tumor  size. 
Presently,  the  reasons  for  decreases 
in  VEGF  expression  in  HuR- 
overexpressing  tumors  are  unclear. 
Our  results  are  highly  reproducible 
when  each  mouse  is  compared  with 
its  cohorts  within  individual  and 
duplicate  experiments.  Furthermore, 
it  is  clone-independent,  in  that  we 
obtained  similar  results  using  two 
different  overexpression  clones. 
However,  we  do  not  yet  know  if  these 
results  are  specific  to  only  ER'  breast 
cancer,  or  also  include  the  ER*  sub- 
type,  as  our  results  deal  primarily 
with  one  cell  line  and  must  be  vali¬ 
dated  in  other  cell  lines  before  we  can 
make  generalized  conclusions. 

It  will  be  interesting  to  test  the 
therapeutic  paradigm  whereby  exog¬ 
enous  HuR  overexpression  in  already 
established  tumors  may  abrogate 
their  development,  perhaps  by  interfering  with  angiogenesis.  If 
correct,  HuR  may  be  involved  in  regulating  a  cluster  of  genes 
involved  in  angiogenesis  and  therefore,  alterations  in  HuR  expres¬ 
sion  may  influence  multiple  downstream  genes  which  operate  in 
blood  vessel  formation.  These  experiments  are  currently  under¬ 
way.  There  has  been  a  great  deal  of  interest  in  blocking  angiogen¬ 
esis,  either  by  increasing  expression  of  anti-angiogenetic  factors, 
such  as  TSP1  or  blocking  pro-angiogenetic  factors  such  as  VEGF. 
It  has  been  challenging  to  work  with  TSP1  due  to  its  large  molec¬ 
ular  size.  The  emerging  evidence  in  cancer  treatment  approaches 
in  this  area,  however,  has  been  that  monotherapies  can  be  more 
readily  overcome  by  cancer  cells,  which  are  constantly  evolving 
in  the  patient.  We  may  have  discovered  a  novel  way  of  limiting 
ER'tumor  growth  by  interfering  with  neo-angiogenesis  at  mul¬ 
tiple  steps  by  modulating  the  expression  of  a  key  RBP  involved 
in  tumor  growth. 

Materials  and  Methods 

Cell  line  and  growth  conditions.  The  MDA-MB-231  cell 
line  was  purchased  from  American  Type  Culture  Collection 
(Manassas,  VA)  and  grown  according  to  the  vendor’s  directions. 

Generation  of  clones  expressing  HA-HuR.  Hemagglutinin 
(HA)  tagged  human  HuR  was  cloned  into  the  Nhel  and 
Xhol  sites  of  the  pZeoSV2  (-)  (Invitrogen™)  vector.  Cells  were 
then  plated  and  transfected  with  either  pZeo  HA-HuR  or 
pZeo  empty  vector  using  Lipofectamine  2000  (Invitrogen™). 
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Figure  S.  Tumors  overexpressing  HA-HuR  have  no  increases  in  apoptosis  but  decreased  blood  vessel  formation  compared  to  control.  (A)  Annexin 
V  staining  reveals  similar  amounts  of  cells  undergoing  apoptosis  between  cells  overexpressing  HA-HuR  and  EV  control  cells.  (B)  Caspase  3  staining 
shows  no  differences  in  the  amount  of  apoptosis  in  the  EV  control  tumors  compared  to  ElA-EluR  tumors  harvested  14  days  post-inoculation.  In  the 
tumors  harvested  on  day  42  post-inoculation,  caspase  3  staining  showed  more  apoptotic  cells  in  the  EV  control  tumors  compared  to  tumors  overex¬ 
pressing  HA-HuR.  CD34  staining  shows  fewer  blood  vessels  in  tumors  overexpressing  ElA-EluR.  (B)  Quantitation  of  blood  vessels  stained  (number  of 
vessels  per  high  power  field  scored)  with  CD34  indicated  significantly  fewer  blood  vessels  in  the  tumors  overexpressing  HA-HuR.  Error  bars  ±  SEM; 
p  <  0.005;  in  photomicrographs  bar  =  27  microns.  Representative  of  n  =  5  sets  of  tumors. 


Selection  was  performed  with  200  |xg/ml  of  Zeocin  anti¬ 
biotic  (Invitrogen™)  and  cells  were  cloned  by  limiting 
dilution. 

SDS-PAGE  and  western  blot  analysis.  Western  analysis  was 
performed  as  described  previously.12  For  detection  of  VEGFa 
and  TSP1  from  tumors,  triple-detergent  RIPA  lysis  buffer  with 
protease  inhibitors  was  used.  Membranes  were  probed  with  anti- 
TSP1  (Abeam)  or  anti-VEGFa  (Abeam)  and  anti- (3 -tubulin 
(Sigma-Aldri+ch).  Secondary  antibodies  used  were:  sheep  anti¬ 
mouse  HRP  (GE  Flealthcare)  or  donkey  anti-rabbit  FIRP  (GE 
Flealthcare)  for  VEGFa  and  TSP1,  respectively.  Proteins  were 
detected  using  chemiluminescence  (GE  Healthcare).  HuR,  TSP1 
and  VEGFa  levels  were  determined  using  Bio-Rad’s  Quantity 
One  software  (Bio-Rad)  normalizing  to  (3-tubulin.  Anti-HuR 
3A2  hybridoma  was  kindly  provided  by  Joan  Steitz  (Yale  School 
of  Medicine). 

In  vitro  growth,  apoptosis  and  cell  cycle  assay.  Fifty  thou¬ 
sand  cells  were  seeded  in  a  24-well  plate.  Cells  were  then  counted 
using  a  hemocytometer  and  trypan  blue  exclusion  dye  24,  48  and 
72  hours  post-seeding  n  =  4  assays.  For  cell  cycle  analysis,  cells 
were  fixed  and  permeabilized,  re-suspended  in  PBS  with  0.2  mg/ 
ml  RNase  A  (Sigma-Aldrich)  and  10  mg/ml  propidium  iodide 
(Sigma-Aldrich),  and  analyzed  on  FACScan  (BD  Biosciences, 
San  Jose,  CA)  and  cell  cycle  analysis  was  performed  using 
Cell  Quest  software  (BD  Biosciences).  Annexin  V  and  7-AAD 
staining  was  performed  using  the  FITC  Annexin  V  Apoptosis 
Detection  Kit  I  (BD  Bioscience)  by  following  the  manufacturer’s 
protocol.  Flow  cytometry  was  performed  on  the  FACScan  (BD 


Biosciences)  and  analysis  was  performed  using  Cell  Quest  soft¬ 
ware  (BD  Biosciences). 

Mice  tumor  inoculations  and  measurements.  Female  athy- 
mic  nude  mice  (aged  8-12  weeks)  were  purchased  from  Harlan. 
For  tumor  inoculations,  100  pL  of  a  1:1  mixture  of  Matrigel 
(BD  Biosciences)  and  RPMI  1640  (GIBCO®)  containing  1 
x  106  MDA-MB-231  cells,  expressing  either  pZeo  HA-HuR, 
pZeo  empty  vector  or  wild-type  clones,  were  injected  into  the 
left  or  right  mammary  fat  pad.  Tumor  volumes  were  calculated 
using  calipers  by  measuring  using  the  formula:  LxWxDx0.5. 
Procedures  were  conducted  according  to  the  University  of 
Missouri  Columbia  Animal  Care  and  Use  Committee. 

Longitudinal  MRI  investigation  and  tumor  volume  analy¬ 
sis.  Magnetic  Resonance  Imaging  (MRI)  was  performed  on  a 
7T/210  mm  Varian  Unity  Inova  MRI  system  (Varian  Inc.,  Palo 
Alto,  CA)  using  previously  published  procedures.44  Mice  were 
anesthetized  with  1-2%  isoflurane  in  oxygen  via  a  nose  cone  over 
the  entire  imaging  period.  Three  mice  were  imaged  weekly  for 
5  weeks  to  monitor  tumor  growth.  Mice  were  imaged  to  obtain 
axial  planes  using  multi-slice  spin-echo  Tl-weighted  (TlW) 
imaging  sequence  applied  with  a  fat-saturation  pulse  to  sup¬ 
press  the  strong  signals  from  fatty  tissues  in  the  chest.  Spin-echo 
diffusion-weighted  imaging  (DWI)  with  b-value  =  1,063  s/mm2 
was  performed  at  week  4  to  assess  the  tumor  tissue  viability,  i.e., 
necrotic  tissue  or  solid  tumor  tissue.  Images  were  obtained  with 
slice  thickness  =  0.8  mm  and  in-plane  resolution  =  0.059  mm  x 
0.156  mm.  The  tumors  were  manually  segmented  using  VnmrJ 
software  (Varian  Inc.,)  to  obtain  the  tumor  volume  in  cm3.  DW 
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images  at  week  4  were  used  to  differentiate  between  necrotic  tis¬ 
sues  and  solid  tumor  tissues. 

RNA  purification  and  real-time  PCR.  RNA  was  extracted 
from  snap  frozen  tumors  in  Trizol  reagent  (Invitrogen™) .  One 
|JLg  of  RNA  was  reverse  transcribed  and  real-time  PCR  was  per¬ 
formed  using  cDNA  (in  triplicate)  for  using  Superscript  III  two- 
step  qRT-PCR  with  SYBR  green  (Invitrogen™).  Primers  were  as 
follows: 

VEGF  sense  5-TTT  CTG  CTG  TCT  TGG  GTG  CAT 
TGG-3'  and  antisense  5 '-ACC  ACT  TCG  TGA  TGA  TTC 
TGC  CCT-3',  TSP1  sense  5-TTC  CGC  CGA TTC  CAG  ATG 
ATT  CCT-3'  and  antisense  5-ACG  AGT  TCT  TTA  CCC  TGA 
TGG  CGT-3',  HIFla  sense  5-TTG  GCA  GCA  ACG  ACA 
CAG  AAA  CTG-3'  and  antisense  5  -TTG  AGT  GCA  GGG 
TCA  GCA  CTA  CTT-3',  GAPDH  sense  5-AGC  CTC  AAG 
ATC  ATC  AGC  AAT  GCC-3'  and  antisense  5-TGT  GGT  CAT 
GAG  TCC  TTC  CAC  GAT-3'.  PCR  reactions  were  performed 
using  Applied  Biosystems  StepOne  PCR  system.  Results  were 
analyzed  using  comparative  CT  method  with  GAPDH  as  endog¬ 
enous  reference  control. 

Microarray.  0.5  |xg  of  total  RNA  was  used  to  make  the  bio¬ 
tin-labeled  antisense  RNA  (aRNA)  target  using  the  Illumina 
TotalPrep  RNA  amplification  kit  (Ambion,  Austin,  TX).  and 
hybridized  to  Illumina  HumanWG-6  (V2_0_R3)  beadchip 
(47,000  genes)  as  per  Illumina  protocols.  After  hybridization,  the 
chips  were  washed  and  stained  with  streptavidin-C3  and  ana¬ 
lyzed  by  BeadArray  reader  (Illumina,  San  Diego,  CA). 

Histology  and  imunohistochemistry.  Tissue  was  rou¬ 
tinely  processed,  formalin-fixed  and  paraffin  embedded  for 
hematoxylin-and-eosin-staining  and  immunohistochemistry. 
Immunostaining  was  performed  using  the  avidin-biotin-per- 
oxidase  complex  method.  Slides  were  incubated  at  room  tem¬ 
perature  with  the  following  antibodies:  anti-cleaved  caspase-3 
antibody  ([2305-PC-100],  Trevigen,  Gaithersburg,  MD); 
anti-CD34  (ab8158,  MEC  14.7,  Abeam,  Cambridge,  MA); 
anti-HIFlct  (NB100-134,  Novus  Biologicals,  Littleton,  CO); 
anti-VEGF  (sc-152  Santa  Cruz  Biotechnology,  Inc.,  Santa  Cruz, 
CA);  and  anti-TSP-1  (clone  A6.1,  MS-420-P1,  Thermo  Fisher 
Scientific,  Fremont,  CA).  Slides  labeled  with  anti-CD34  or 
TSP-1  were  incubated  with  a  biotinylated  secondary  antibody 
DAKO  for  CD-34  and  rabbit  anti-mouse  IgG  (for  TSP-I)  fol¬ 
lowed  by  streptavidin-linked  HRP  (DAKO).  Cleaved  cas- 
pase-3  and  VEGF  slides  were  incubated  with  a  HRP  antibody 
(EnVision™,  DAKO).  Bound  antibodies  were  visualized  with 
peroxidase  substrates:  DAB  3,  3'-diaminobenzidine  solution 
(DAKO)  or  NovaRED™  (Vector  Labs,  Burlingame,  CA).  For 
the  TUNEL  apoptosis  detection  assay  we  used  the  In  Situ  Cell 
Death  Detection  Kit,  TMR  red  (Roche  Applied  Science)  and 
followed  the  manufacturer’s  protocol  for  paraffin  sections. 

Statistics.  All  error  bars  represent  standard  error  of  the  mean. 
The  p  values  were  calculated  using  the  two-tailed  Student  t  test. 

Statistical  analysis  of  microarray  data.  Analysis  of  microarray 
gene  expression  data  was  primarily  performed  using  the  Linear 
Models  for  Microarray  Data  (limma)  package45  and  the  lumi 


package,46  available  through  the  Bioconductor  project47  for  use 
with  R  statistical  software.  Quantile  normalization  was  used  for 
between  chip  normalization.  Statistical  analysis  was  performed 
usingmoderated  t-statistics  applied  to  thelog-transformed(base2) 
normalized  intensity  for  each  gene.  Because  two  measure¬ 
ments  were  taken  from  each  mouse,  the  dependency  between 
paired  measurements  was  accounted  for  by  a  modified  mixed 
linear  model  that  treated  each  animal  as  a  block.  The  within- 
block  correlations  were  constrained  to  be  equal  between  genes, 
and  then  information  was  borrowed  across  genes  to  moder¬ 
ate  the  standard  deviations  between  genes  via  an  empirical 
Bayes  method.  The  contrast  of  interest  computed  and  tested 
was  the  difference  between  overexpressor  and  control  vector. 
Adjustment  for  multiple  testing  was  made  using  the  false  dis¬ 
covery  rate  (FDR)  method  of  Benjamini  and  Hochberg.48  We 
chose  10%  as  our  FDR-cutoff  for  declaring  statistical  signifi¬ 
cance,  and  a  threshold  of  at  least  three-fold  (up  or  down)  for 
declaring  a  biologically  significant  change  in  expression.  Gene 
ontology  (GO)  analyses  were  carried  to  test  the  association 
between  Gene  Ontology  Consortium  terms  and  the  list  of  dif¬ 
ferentially  expressed  genes.  In  defining  the  gene  universe  for 
the  analysis,  non-specific  filtering  was  used  to  increase  statis¬ 
tical  power  without  biasing  the  results.  This  filtering  selected 
only  probes  on  the  Illumina  array  which  had  both  an  Entrez 
gene  identifier,49  a  GO  annotation  (as  provided  in  the  lumiHu- 
manAll.db50  annotation  data  package  and  GO.db51  annotation 
maps)  and  an  interquartile  range  of  >0.1  on  log2  scale  across  all 
samples.  Using  this  gene  universe,  GOstats52  was  used  to  carry 
out  conditional  hypergeometric  tests  which  exploit  the  hierar¬ 
chical  nature  of  the  relationships  among  the  GO  terms  for  con¬ 
ditioning.53  We  carried  out  GO  analyses  for  over-representation 
of  biological  process  (BP),  molecular  function  (MF)  and  cel¬ 
lular  component  (CC)  ontologies  and  computed  the  nominal 
hypergeometric  probability  for  each  GO  category.  These  results 
were  used  to  assess  whether  the  number  of  selected  genes  associ¬ 
ated  with  a  given  term  was  larger  than  expected  and  a  p-value 
cutoff  of  1%  was  used.  GO  categories  containing  less  than  10 
genes  from  our  gene  universe  were  not  considered  to  be  reliable 
indicators  and  are  not  reported. 

RNA  immunoprecipitation.  RNA  immunoprecipitation 
was  performed  as  previously  described.42 
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